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The absorption spectrum of phenol as observed with a 3-m grating spectrograph consisted 
of 330 bands. A tentative analysis is made on the basis of the fundamental ground-state 
frequencies 245, 397, 527, 811, 824, and 999 and the fundamental excited-state frequencies 
208, 374, 476, 562, 783, 935, 975, 1274, and 1565. These agree well with Raman data. 





INTRODUCTION 


HE substitution of the OH group on the 
benzene ring modifies the electronic struc- 
ture of the ring by an inductive effect and by a 
resonance effect involving the unshared pairs of 
electrons on the oxygen atom. This results in an 
intensification of the spectrum and an increase 
in the number of allowed vibrational transitions 
over that of benzene similar to the situation 
found in aniline! and chlorobenzene.? The ab- 
sorbing power of phenol has been estimated to 
be 230 times that of benzene.* 

The ultraviolet absorption spectrum of phenol 
has been investigated by Savard,’ Masaki,® 
Asagoe, Shimokawa, and Kageyama,* and by 
Hall.* These investigators have not presented the 
complete spectrum nor have they presumably 
resolved the doublet fine structure of phenol. 


1N. Ginsburg and F. A. Matsen, J. Chem. Phys. 13, 
167 (1945). 
a9 H. Sponer and S. H. Wollman, J. Chem. Phys. 9, 816 

41). 

3K. Asagoe, Y. Shimokawa, and S. Kageyama, Jap. J. 
Phys. 14, 11 (1940). 

4J. Savard, Ann. de Chim. 11, 287 (1929). 

5K. Masaki, Bull. Chem. Soc. Japan 11, 346 (1936). 

6M. B. Hall, Phys. Rev. 61, 736 (A) (1942). 


Savard analyzed 120 bands on the basis of four 
frequencies. However, he assumed the vibra- 
tional frequencies to be the same in the excited 
state as in the ground state which invalidated his 
assignment, and further, he did not make corre- 
lation with Raman data: Neither Masaki nor 
Hall give any detailed analysis. Asagoe, Shimo- 
kawa, and Kageyama made an analysis on the 
basis of five fundamental frequencies in the 
excited state and two in the ground state. 

The Raman spectrum has been studied by 
Kahovec and Reitz,’? Masaki,®> and Kohlrausch 
and Wittek.® 


EXPERIMENTAL 


The absorption spectrum of phenol was studied 
in the first order of a spectrograph of the Eagle 
type, with a 3-meter grating ruled 15,000 lines 
per inch, to determine the optimum conditions. 
A series of plates were taken in the second order 
at temperatures of the absorption call ranging 
from 24°C to 120°C while the temperature of the 


7L. Kahovec and A. W. Reitz, Monatsh. 69, 363 (1936). 
8K. W. F. Kohlrausch and H. Wittek, Monatsh. 74, 1 
(1943). 
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Fic. 1. Pressure development of the phenol spectrum. Cell temperature 24°C; reservoir temperatures: 
(A) —37°C, (B) 0°C, and (C) 24°C. 


phenol reservoir, which controlled the pressure, 
was varied between —37°C and 90°C. The ac- 
companying figure (Fig. 1) illustrates the devel- 
opment of the spectrum. 

The plates used were Eastman Kodak 103-0. 
The phenol, “Special Absolute Pure’’ Phenol 
C. P., prepared by the National Aniline and 
Chemical Company, was redistilled under vac- 
uum. The absorption cell, as all quartz, was 70 
cm long. The light source was a 500 volt 2.5 kva 
hydrogen discharge tube, sufficiently powerful to 
give exposures in two hours even with the use of 
the necessary filters. 

No photo-decomposition was observed for 
phenol. The absorption cell was pumped out and 
refilled several times during the course of the 
investigation, but there was no_ noticeable 
decrease in the intensity of the exposure durings 
several runs. 

Over 300 bands were observed in the spectrum 
in the region 2490-2830A. Table I gives the 
wave-length, wave number, intensity, and as- 
signment of these bands. Since many of the bands 
appear as very sharp doublets, measurements 
were made to the centers of absorption in these 
cases, and the method was carried over to the 
broad, more diffuse bands. For the doublets, the 
precision was estimated to be 0.3 cm=, with the 
accuracy diminishing to 2 cm~ for the weak, 
diffuse bands. In the table, all doublets are 
grouped in pairs. In making intensity assign- 
ments, allowance was made for the effect of 
varying pressures. 


In the spectrum of phenol some seventy 
doublets were found with a mean separation of 
three wave numbers. This doublet character}for 
a number of bands was found to depend on the 
temperature and pressure of the phenol. Some 
doublets were found only at high temperatures 
and pressures; others, which were broad and 
diffuse at both low and high temperatures and 
pressures, were doublets at intermediate values. 
Almost without exception, however, the iden- 
tifiable fundamental frequencies produced per- 
sistent doublet bands. 





Fic. 2. Enlargement (X12) of O—O band in phenol 
showing doublet structure. 
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TABLE I. The spectrum of phenol. The intensities are designated as vw =very weak, w=weak, m=median, ms =medium 
strong, s=strong, vs=very strong, and ws =extremely strong, B=broad and D=diffuse. (Bands grouped 
together are the components of a doublet.) 














Wave- Fre- Wave- Fre- 
length quency Inten- length quency Inten- 
(A) (cm~) sity Assignment (A) (cm~) sity Assignment 
2827.93 35,351.2 vw 0—999 2766.64 36,134.3 wwBD 0-—999+783 
19.09 35,462.0 vw 0—824-—61 65.89 36,144.0 ww 0—245+40 
18.76  35,466.2 65.62 36,147.6 
65.37 36,150.9 
14.11 35,5248 w 0—824 
13.83 35,528.3 65.05 36,155.0 ww 
64.86  36,157.5 
13.04 35,538.33 wwB 0-811 64.63 36,160.5 0—397+208 
07.78 35,604.9 vw 0—527—221 63.27 36,178.2 ww 


63.08  36,180.8 


2800.55 35,696.8 wwB 
62.61 36,187.00 ww 





2798.28 35,725.7 wwB 0-—999+374 62.42 36,189.5 
nty 95.62 35,759.5 ww 0—527-61 62.11 36,193.5 w 0—527+374 
1 of 95.34 35,763.3 61.89  36,196.4 
“ 
,for 94.68 35,771.8 ww 61.28 36,2044 w 
the 94.43 35,774.9 61.09  36,206.9 
ome 93.63 35,785.2° wwB 60.34 36,216.7 vw 
ires 
aver" 92.46 35,800.2 wwD 60.10 36,219.9 ww 
and 91.83 35,808.3 wwB 59.78 36,2241 w 0-2x61 
; 59.57 36,226.8 
ues. 91.15 35,8170 w 
len- 58.63 36,239.2 vwwB 
90.78 35,822.5 m 0-527 0—999+476 
per- 90.50  35,825.3 58.18  36,245.1 ww 
85.27 35,892.6 ow 0—397-61 57.96 36,248.0 vwB 
83.73 35,912.5 vow 0—999+562 0—811+374 57.59 36,252.9 ms 
83.54 35,914.9 57.38  36,255.6 
80.71 35,951.5 m 0—397 55.50 36,280.4 ow 
55.27 36,283.4 
79.00 35,973.6 vw 
54.95 36,287.6 ws 0-61 (0—999+935) 
76.08 36,011.4 ow 0—811+476 54.74  36,290.7 
75.89  36,013.6 
54.12 36,298.5 m 0—527+476 
75.23 36,022.4 vw 53.91  36,301.3 
74.76 36,028.5 ww 0—527+208 53.18 36,3109 s 0—824+783 0—245+208 


52.96  36,313.8 
73.92 36,039.4 vw 
73.71 36,042.2 50.82 36,342.1 vs 
50.62  36,344.7 


72.66 36,055.88 vw 
50.29 36,348.6 wus 0-0 





71.74 36,067.8 ww 0—221—61 50.08  36,351.9 
71.53 36,070.5 
49.30 36,362.2 vw 
70.19 36,088.0 ww 0—824+ 562 49.12 36,365.1 
69.99  36,090.6 
47.72 36,383.1 ww 
68.95 36,105.4 ww 0—245 
68.74  36,106.9 47.54 36,385.5 ww 0—527+562 
67.14 36,127.7 s 0-221 47.27 36,389.00 s 0+40 


| 66.94 36,130.6 47.06 36,391.8 
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TABLE I.—Continued 








Wave- Fre- Wave- Fre- 





length quency Inten- length quency Inten- 
(A) (cm~) sity Assignment (A) (cm~!) sity Assignment 
2746.59 36,398.00 m 2716.19 36,805.4 wwB 0+2xX208+40 
46.39  36,400.7 


15.53 36,814.3 wD 0—811+1274 
15.23 36,8184 wD 


45.39 36,415.1 vw 
44.48 36,4260 wwD 0-397+476 


14.75 36,824.9 vs 0+476 
43.38 36,440.6 vw 14.58 36,827.2 
41.77 36,462.0 wwB 0-—824+935 12.72 36,852.55 wwB 0-—61+562 
41.04 36,471.7 vw 11.85 36,864.3 vw 0+476+40 
11.68  36,866.6 
39.27 36,495.33 wBD 0—-61+208 
11.26 36,872.3 ww 
38.68  36,503.2 vw 0—824+975 0—221+374 11.08  36,874.7 
38.47 36,506.0 
10.12 36,887.8 vw 0—245+783 0—397+935 
36.99  36,525.7 vw 09.94  36,890.3 
36.77 36,528.6 
; 08.39 36,9114 m 0+562 0—221+783 
36.05 36,538.3 ww 08.20 36,914.0 
35.84  36,541.0 


07.14 36,9284 vw 0—347+975 0+208+374 
06.00 36,944.0 vw 


35.08  36,551.2 vw 
34.66 36,556.8 ms 0+208 


34.47 36,559.4 04.80 36,960.4 wwB 
32.73 36,582.6 wwB 0—245+476 04.09 36,969.8 wwB 0+3X208 
31.50 36,599.1 wwB 0+208+40 03.50 36,978.1 vwB 
31.01 36,605.7 wD 0—221+476 0—527+783 02.87 36,985.9 ow 
30.79  36,608.8 02.72  36,988.8 
29.58  36,624.9 wwB 0-—999+1274 01.42 37,006.66 ww ° 
01.22  37,009.4 
: 631.2 
_—— es wee 2699.34 37,035.1 w 0-245+4.935, 0+208-+476 
wal _ 0—61+2X374 
26.73 36,663 i vw 0—245+562 0—61+374 99.10 37,0384 
: ee 98.03 37,053.1 vw 
25.89 36,6744 w : a 
253.71 36,676.9 97.06 37,066.4 wD 0—221+935 
96.75 37,070.7_ s 0—61+783 
23.66 36,704.5 vow 96.58 37,073.0 
222.92 36,714.3 vow 03 37,080.6 0—245+975 
22.63 36,718.3 ™ a ‘os + 
95.53 37,087.5 wwB  O0—824+1565 
22.31 36,722.7. m 0+374 
22.09  36,725.6 94.99 37,094.9 m 0+2X374 O—527+1274 
94.77 37,097.9 
20.52 36,746.7 ww 
93.28 37,1185 wD 0+208+562 
19.73 36,757.5 wwB 0—527—935 0—811+1565 
19.25 36,764.0 ww 0—61+476 92.75 37,125.7 msD 
19.07 36,766.4 0+374+40 0+2%X208 92.59 37,128.0 
18.42 36,775.2 ww 92.27 37,131.5 ws 0+783 0+2X374+40 


92.10 37,134.7 


17.51 36,787.5 vow 
89.50 37,170.66 m 0+783+40 
17.09 36,793.2 vw 0—527+975 89.29 37,173.5 


















935 


374 


476 


74 
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TABLE I.—Continued 


PHENOL VAPOR 








Wave- 


Fre- 


Wave- 


Fre- 








length uency Inten- length quency Inten- 
(A) cm?) sity Assignment (A) (cm~) sity Assignment 
2688.62 37,182.8 wwD 2667.65 37,475.0 ww 0+2x 562 
87.89  37,192.9 w 66.18 37,495.7 wB  0+4+208+374+562 
87.67 37,195.9 0+374+476 = 042081935" 
86.20 37,216.3 m 0+783+240 
65.36 37,507.2 ww 043744783 
86.03  37,218.6 65.18  37,509.7 0+208+2X476 
85.75  37,222.5 0-614+935 0—397+1274 
85,50 372260 - . 64.54 37,5188 ww  0—-397+1565 
84.86 36,2348 w 63.90 37,527.8 wwB  0+208+975 
84.62  37,238.2 
as meee o 62.83 37,542.9 ww 0—614+476+783 
83.67 37,251.3 61.57 37,560.7 ww 0—61+41274 
82.80 37,2634 mBD 0—61+975 61.37 37,563.5 
$2.22 37,271.5 mBD 60.60 37,574.3 vwwB 0+2X374+476 
0—527+476+1274 
81.91 37,275.88 wD 0—61-+208+783 
59.83 37,585.2 vwD 
81.30 37,284.3 ws 0+935 043744562 
81.11  37.286.9 59.17 37,594.5 ww 0—61+3744935 
80.90 37,289.8 0+208+476+562 
80.51  37,295.3 ms 57.87 37,613.0 w 0+476+783 
80.30 37,298.2 57.65  37,616.0 
79.65  37,307.2 m 0+2X374+208 57.18 37,622.7. s 0+1274 
79.48  37,309.6 57.01 37,625.1 
79.16 37,314.1 wBD 56.24 37,636.0 ww 0—61+562+783 
0—614+374+975 
78.47 37,323.7 s 0+975 0+935+40 
78.22  37,327.2 0+2%X 2084562 54.51 37,660.5 ww 0+1274440 0+374+935 
0—61-+476+542 54.37 37,662.5 0+2% 3744562 
1742  37,338.3 w 0+208-+783 52.06 37,695.3 wmwBD 0—221+1565 
77.27 37,340.4 0+562+783 0+374+975 
76.79 37,347.11 ww 47.42 37,7614 ms 0+476+4+935 
47.26 37.763.7 0+3744476+562 
76.15  37,356.0 ww 
46.90 37,768.8 ww 0—61+208+1274 
75.51 37,365.0 ww 0+975+40 
45.95 37,782.4 ww 0—61+4+562+935 
74.63 37,377.3 wwB 0—245+1274 
45.54 37,788.2 ww 
73.77 37,389.3 vwB 0—527+1565 
0+476+562 44.80 37,798.8 vow 0+476+975 
73.39 37,394.6 vow 44.01 37,810.1 ww 
72.81 37,402.7 wwB  0—221+1274 43.31 37,820.1 ww  0+208+476+783 
70.62 37,433.4 wB  0—61+208+935 42.71 37,828.7 ww  0+208+1274 
69.59 37,447.88 ww 0—61+374+783 41.64 37,844.0 ww 045624935 
0+374+2X562 
68.94 37,456.9 www 
68.78 37.459.2 41.17 37,850.7. w 0—61+1565 
40.97 37.853.6 
68.25 37,466.6 vwwD 
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TABLE I.—Continued 
Wave- Fre- Wave- Fre- 
length quency Inten- length quency Inten- 
(A) (cm~) sity Assignment (A) (cm~) sity Assignment 
2640.50 37,860.3 vw 0+2X476+362 2615.71 38,219.1 msD 042x935 
15.55  38,221.5 
39.31 37,877.4 vwBD 
14.80 38,2324 wwB 0424764935 
38.87 37,883.7 ww 0+562+975 
; 0+2X374+783 13.84 38,246.5 vwBD 
38.40 37,890.5 ww 12.90 38,260.2 wD  0+4+2x935+40 
0-+935+975 
37.83 37,898.5 wwB  0+208+562+783 
11.76 38,277.00 wwB 0+2X476+975 
36.74 37,914.3 s 0+1565 0+2X783 0+208-+783+935 
36.64 37,915.8 
10.70 38,292.5 wwB 0437442783 
36.12 37,923.3 vw 
09.66 38,307.7 vwwB 0+208+476+1274 
34.50 37,946.6 wwBD 0+476+2X562 
08.42 38,326.0 wwB  04476+562+935 
33.44 37,961.6 ww 
07.20 38,343.9 ww 0+61+783+1274 
32.73 37,972.1 ww 0+208+476+935 
05.58  38,367.7 ww 0+2X374+1274 
30.99  37,997.2 ww 0+374+41274 
03.57 38,397.3 wBD  0+4208+562+1274 
30.48 38,004.6 »vBD 0+208+476+965 
02.98 38,4060 msB  0+783+1274 
29.70 38,015.8 ww 0+2X 562+935 
28.65 38,031.0 ww 0+2X374+935 00.40 38,444.1 vow 0+374+783+935 
27.85 38,042.6 ww 0—61+783+935 Pray ore 
26.71 38,0591 wD  0—61+208+1565 — ee 
26.32 38,064.7  msD 93.57 38,545.4 ww 0+476+783+935 
26.12  38,067.7 vs 0+783-+935 93.21 38,550.7 
25.86 38,071.4 92.76 38,557.4 m 0+935+1274 
25.49 38,0768 w 92.54  38,560.7 
= ee 92.04 38,568.1 ww 
24.65 38,089.0 ww 
90.39 38,592.7 ww 0437442935 
24.24 38,094.9 ww 0+208+562+975 90.18  38,595.8 0+975+1274 
23.80 38,101.3 vw 0+476-+1274 88.12 38,626.5 ww 0+562+783-+935 
23.51 38,105.5 m 0+783+975 87.66 38,633.4 ww 0+3744+935+975 
23.37 38,107.4 
83.55 38,695.9 wBD 043X783 
22.55 38,119.6 wwB 0+208+1565 0+476+2X935 
0+208-+2X783 
81.24 38,729.5 wwD 
21.68 38,132.1 wwB 0437444764935 
78.38 38,772.4 wwD 0+476+2x975 
20.75 38,145.7 vow 
77.36 38,787.8 vwBD 
19.99 38,156.7 wwB 0—61+2X935 
73.29 38,849.1 wBD 04935+1565 
19.10 38,169.7 wwB  0+4+476+562+783 
72.74 38,857.4 vow 
18.16 38,183.4 vow 0+562+1274 
72.20 38,865.6 ww 0+2X476+1565 
17.17 38,197.8. ww 0—61+935+975 
71.58  38,875.0 vw 
38,213.7 
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TABLE I.—Continued 








Fre- Wave- Fre- 
quency : length quency , 
(cm~) Assignment (A) (cm~}) Assignment 


38,886.6 0+975+1565 2539.54  39,365.4 0+476+-975+1565 


38,894.7 0+2X 1274 38.85  39,376.1 0+476+2X 1274 
0+783+975+ 1274 








38,909.5 
34.44 39,444.6 


32.10  39,481.0 0+4X783 0+2X1565 
39,000.5 0+783+2X935 
22.67  39,628.6 0+762+935+ 1565 
20.14 39,668.4 0+783+-975+ 1565 
39,040.0 0+2X 562+ 1565 bea didi 
0+783+935+975 13.02 39,780.8 0+2X935+1565 
39,151.6 0+3X935 10.42 39,822.0- 0+935+975+ 1565 
39,187.0 0+1274+1565 03.30 39,935.2 
39,329.2 01.01 40,080.7 


eee 0+783+935+1274 2491.65 40,122.0 0+935+1274+1565 
9,342.2 


38,936.7 


39,026.0 0+476+935+1274 











TABLE II. Fundamental vibrational frequencies of phenol. 








Ground state Excited state 
Ultraviolet Raman 


Asagoe, Asagoe, 
Shimokawa, This Kahovec, Kohlrausch, This Shimokawa, 
Kageyama research Reitz Wittek research Kageyama 


(235) 245 w 242 (6b) 240 (5) 208 ms 

397 m 374 m 
(526) 527 m 532 (5) 532 (5) 476 wus 473 

(614) 623 vow 617 (5) 617 (8) 562 m 
811 wwb 812 (8b) 812 (100) 783 vvs 779 
824 w 826 (3) 829 (12) 
Reported as 999 ww 1000 (2) 1001 (10) 935 vvs 932 
differences 1025 (8) 1026 (7) 975 vus 982 
1253 (4b) 1250 (26) 1273s 1273 

1598 (5d) 1595 (6) 1565 s 














DISCUSSION 999 cm! were found corresponding to the sym- 
metrical Raman frequencies 242, 532, 812, and 
1000 cm reported by Kahovec and Reitz 
(Table II). In addition, a band of medium inten- 
sity was found 824 cm to the red of the 0—0 
band which probably corresponds to the anti- 
symmetric Raman frequency 826 cm™. If this 
association is correct, it is the first instance in 
our work of an antisymmetric vibration entering 
into the spectrum in this way. The 617 cm™ 
Raman frequency did not appear as 0—1X»; 

® y will be used to represent fundamental frequencies in transition although Masaki reported a band 604 
the ground state, and »’ in the excited state. cm~! to the red from 0—0, and Asagoe, Shimo- 


The 0—0 transition, the most intense band in 
the spectrum, appeared as a doublet (see Fig. 2), 
the wave numbers of the components being 
36,348.6 and 36,351.9 to be compared with the 
frequencies given for a single band of 36,351 cm 
reported by Hall and 36,356 reported by both 
Masaki, and by Asagoe, Shimokawa, and 
Kageyama. 

The 0—1 Xv»; ° transitions 245, 526, 811, and 
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kawa, and Kageyama reported 614 as an im- 
portant difference frequency. 

The medium band 397 cm to the red from 
0—0 has also been taken as 0—1 Xv even though 
it is missing in the Raman spectrum. Other 
monosubstituted benzenes have a strong Raman 
line in this neighborhood, and in fact, Kahovec 
and Reitz leave a space for this frequency in their 
table. In addition, this frequency combined well 
with a number of frequencies in the excited state. 
It was thought to be the same mode of vibration 
as the strong frequency 374 of the excited state. 

The low intensity of 623, 811, and 999 cm 
even at high phenol pressures and high cell tem- 
peratures seemed to indicate some selection rule 
against their appearance. All three have a greater 
Raman intensity (8, 10, and 10, respectively) 
than 532 (intensity 5) which appeared strongly 
even at low phenol pressures and low cell tem- 
peratures. That the Boltzmann factor does not 
account for this low intensity may be seen from 
the ratios of the Boltzmann factors of the three 
higher frequencies to the Boltzmann factor of the 
532 frequency; these are 0.67, 0.26, and 0.11, 
respectively, at a temperature of 27°C. 

In addition. to the frequencies which may be 
correlated with the Raman data there appeared 
two strong bands 61 and 221 cm to the red 
from the 0—0 band. These may be difference 
frequencies and because of the high intensity 
probably involve a low frequency in the ground 
state. These combined with a number of fre- 
quencies in the excited state. Hall reported 
values of 60 and 218 cm while Asagoe, Shimo- 
kawa, and Kageyama reported 60 and 220 cm“. 
~ The identification of fundamental frequencies 
in the excited state was difficult because only an 
indirect correlation with the Raman data was 
possible and because of the large number of bands 
in the phenol spectrum. The bands were so closely 
spaced that if frequencies were chosen at random, 
the probability of coincidence (3 cm™@') was 
between 0.4 and 0.5. 

The fundamental frequencies in the excited 
state were chosen on the following basis: the 
intensity of the 0+1X >,’ transition, the correla- 
tion with Raman and ground-state frequencies, 
and the combinability with other chosen fre- 
quencies. These frequencies are listed in Table II 
opposite the Raman frequencies with which they 
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are probably associated. The association was 
made with all the strong totally symmetric fre- 
quencies with the exception of 1274. Asagoe, 
Shimokawa, and Kageyama have associated the 
1272 excited state frequency (1274 in this 
research) with the Raman line, 1252. This is not 
an impossible association, but it is generally 
believed that electronic excitation loosens the 
molecule so that the vibration frequencies of the 
excited state would be smaller than those of the 
ground state. 

A strong band appeared 40 cm to the blue 
from the 0-0 band. This combined strongly with 
all the fundamental frequencies in the excited 
state and with many combinations of these fun- 
damental frequencies. It did not, however, com- 
bine withany of the frequencies in the ground state. 

Due to the fact that the assignment of funda- 
mental frequencies (either the choice of the 
frequency or its precise value to a single wave 
number) was not absolutely certain, that all 
possible combinations of the fundamental fre- 
quencies did not appear, and that a high prob- 
ability for accidental coincidences existed, the 
situation with regard to the combination bands 
was unsatisfactory. In this research a large 
number of observed frequencies were left unas- 
signed because the authors could not be certain 
from the data that transitions of the type 
0-—1Xv:+1Xvj/+1X»;’ were present. Some as- 
signments of the type 0+1X»v/+1Xv/+1X»:’ 
were included even though a number of this type 
were missing. 

All two-quantum number transitions, 0+2 X »;’ 
and 0+1X»v,’+1Xv; were present with the 
exception of 0+526+40. A large number of the 
three-quantum number transitions were missing. 
Among the better series were 0+783+935+ 
1Xv;' of which no members were missing, and 
0+2X935+1 Xv,’ which did not contain 208, 
476, and 562. 

It would seem that further work on the mono- 
substituted benzenes should take away some of 
the uncertainty in making the assignments and 
lead to the formulation of definite selection rules 
if they exist. Research of such character is in 
progress in this laboratory. 

The authors acknowledge with gratitude, the 
continued support of the University of Texas 
Research Institute. 
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The Structure of Liquid Argon 


Louis H. Lunp 
University of Missouri, Columbia, Missouri* 


The method of Wall for calculating the free volume per atom in a liquid has been applied to 
liquid argon, using the atomic distribution curves at 84.3°K and at 91.8°K given by Eisenstein 
and Gingrich. From the free volumes at these temperatures, the entropy of the liquid at the 
boiling point is calculated to be 15.8 cal. degree! mole~!. This value, together with the cal- 
culated entropy, 32.2 cal. degree~! mole, of argon vapor at the boiling point, give an entropy 
of vaporization of 16.4 cal. degree! mole“. The latent heat pf vaporization is thus calculated 
to be 6.01 kilojoules per mole; this is to be compared with the observed value of 6.3 kilojoules 
per mole. Presented also is a comparison of theoretical and experimental atomic distribution 
curves at the temperatures, 84.3°K and 91.8°K. 





1. INTRODUCTION 


NE of the most fundamental characteristics 
of a liquid is its structure. Compared to the 
structure of a crystal, the structure of a liquid 
is much less evident, much less exact over ex- 
tended regions, and in a microscopic sense much 
less simply specified. Perhaps the best description 
of the structure of a liquid that can be given is a 
“statistical’’ one in terms of its atomic, or 
molecular, distribution curve! which gives the 
average number of atoms between distances R 
and R+dR from any arbitrarily chosen atom in 
the body of the liquid. The ultimate utility of 
these curves arises from the fact that they present 
a graphical representation of time-averaged 
atomic or molecular arrangements induced under 
the action of characteristic interatomic or inter- 
molecular forces. 

Several different approaches have been at- 
tempted to determine these distribution curves, 
and some attempts have been made to use them 
in deducing macroscopic properties of elementary 
substances. They have been determined? directly 
from experimental data, and they have been 
obtained* from several theoretical approaches 
which rely in different ways on experimental data 
for the evaluation of otherwise unknown param- 


* Now with the Laboratory, Lukas-Harold Corporation, 
operators of the U. S. Naval Ordnance Plant, Indianapolis, 
Indiana. 

1 P. Debye and H. Menke, Physik. Zeits. 31, 379 (1930). 

2 N.S. Gingrich, Rev. Mod. Phys. 15, 90 (1943). 

3 J. G. Kirkwood and E. M. Boggs, J. Chem. Phys. 10, 
394 (1942); J. Corner and J. E. Lennard-Jones, Proc. Roy. 
Soc. A178, 401 (1941); G. S. Rushbrooke and C. A. Coul- 
son, Phys. Rev. 56, 1216 (1939); C. N. Wall, Phys. Rev. 
tea (1938); J. A. Prins and Petersen, Physica 3, 147 
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eters. Thus far, the diffraction of x-rays by liquid 
elements has supplied? the majority of known 
distribution curves for relatively simple liquids, 
and some of these* have been used in simplified 
calculations to obtain quantities such as heats of 
fusion or vaporization. 

The liquids sodium, potassium, and argon have 
been the subjects of special studies. Kirkwood 
and Boggs® applied a general method to obtain 
an approximate solution for the radial distribu- 
tion function of a system of ‘‘hard spheres ;” they 
then compared this to experimentally obtained 
distribution curves for argon. Rushbrooke® cal- 
culated the shape of the first peak on the dis- 
tribution curve for liquid argon on the basis of 
the Lennard-Jones and Devonshire theory’ of the 
liquid state of inert gases. In the present work, 
Wall’s theory* has been applied to curves ob- 
tained by Eisenstein and Gingrich® for liquid 
argon at T=84.3°K and T=91.8°K. 


2. THEORY 


Wall’s theory® assumes a simple model of the 
quasi-solid type in which each atom is trapped 
by its immediate neighbors in a small spherical 
cell. The potential energy function within the 
cell is assumed to be a discontinuous ‘‘box’’ type 


4*C. N. Wall, Phys. Rev. 54, 1062 (1938); N. S. Gingrich 
and C. N. Wall, Phys. Rev. 56, 336 (1939); J. H. Hilde- 
brand, J. Chem. Phys. 7, 1 (1939). 

5 J]. G. Kirkwood and E. M. Boggs, J. Chem. Phys. 10, 
394 (1942). 
6G. S. Rushbrooke, Proc. Roy. Soc. Edinburgh 60, 182 

1940). 

7 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A163, 53 (1937). 

8C. N. Wall, Phys. Rev. 54, 1062 (1938). 

°A. S. Eisenstein and N. S, Gingrich, Phys. Rev. 62, 
261 (1942), 
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Fic. 1. A comparison of theoretical and experimental 
atomic distribution curves for 4rRp(R) vs. R for liqud 
argon at 84.3°K. Solid line is experimental. Dotted line is 
theoretical (¢ by R—C method). Circles are theoretical 
values (¢ by iteration). 


in which the potential is practically constant 
throughout the cell except at the cell boundary 
where it becomes infinite. This implies a uniform 
probability of finding the atom within equal 
volumes of the cell. The analytical expression he 
derives for the distribution function is 








3NiR 5(Ri—R)? 
4rR*p,(R) = | 1- 
oR, 4c? 
|Ri-R| |Ri—R|? 
x} 1-— + (1) 
20 400° 
for 
R,-20=RE=R i426, 
and 


42 R’pi(R) =0 


for all other values of R. 

In this expression, pi(R) is the density of 
nearest neighbors surrounding an atom in the 
liquid, i.e., the density of those atoms whose cell 
centers lie on the first co-ordination sphere of 
radius Ri, 42R?p;(R)dR is the number of nearest 
neighbors between distances R and R+dR from 


a given atom, MN; is the average number of 
nearest neighbors of a given atom, ¢ is the radius 
of the sphere whose volume represents the free 
volume for the trapped atom (i.e., its center) 
within it. It is assumed that the density of 
nearest neighbors, pi(R), is practically the total 
density of atoms in the interval O=R=R;, and 
hence that in this interval, Eq. (1) is the atomic 
distribution function for the liquid. 

If it is desired to extend the range of the dis- 
tribution curve, it is necessary to consider 
density contributions from atoms whose cell 
centers lie on co-ordination spheres of order 
higher than the first. If p; is the density con- 
tribution of the V; atoms lying on the co-ordina- 
tion sphere of radius R;, then the total distribu- 
tion curve is given by 

4nR?’p(R)=4rR? D p:(R), (2) 
i=l 
where p;(R) is given by an equation of the same 
form as (1) with NV, and R; replaced by N; and R,. 
For the present, we shall restrict the calculations 
to the interval O=R=R3. 


3. APPLICATION TO LIQUID ARGON 


The parameters o, Ni, Ri are found by fitting 
the distribution function (1) to the experimental 
values for liquid argon. Using for convenience, 
4rRpi(R) rather than 47R?pi(R), Eq. (1) may be 
written in the form 


yi  5(Ri—R)? 
x 





= 


wm-y 4 


|Ri—R| |Ri-—R|? 
{1 + | (3) 
20 40a? 


where y=427Rp,(R) and y:=3Ni/5eR; and is the 





TABLE I. Calculated values of o by iteration. 














T =84.3°K T =91.8°K 
R (A) 47Rpi(R)  o@ (A) R (A) 42Rpi(R)_ o (A) 
: exp (calc) (exp) (calc) 
3.4 1.07 0.34 3.2 0.37 0.43 
’ jon 2.09 0.37 3.3 0.82 0.44 
3.6 2.77 0.36 3.4 1.31 0.44 
3.7 3.31 0.42 
Average 0.37A Average 0.44A 
R,=3.75A R,=3.74A 
y1=3.36 41 =2.48 
N,=7.8 atoms N,=6.8 atoms 
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maximum value of y. From the experimental 
values of 42R?’pi(R), a table of y vs. R may be 
constructed and y; and R; are thus immediately 
found. « may then be found by solving (3) by 
iteration for several paired experimental values 
of Rand y. Table I shows the experimental values 
of 4rRpi(R) and the calculated values of o for 
several values of R<R, for liquid argon at 
84.3°K and 91.8°K. Clearly the agreement 
between the values of o calculated in this 
manner at each temperature is quite good. 

Rushbrooke and Coulson!® (hereafter referred 
to as R—C) have proposed two somewhat sim- 
plifying relations to determine the parameters o 
and N;. They are 


Ni/Ni=Riyi/ Ri, (4) 
1+Ni+(N2/2) = (4/3)aRi' po, (S) 


where y, is the value of 4r%Rp(R), corresponding 
to R;, and po is the average macroscopic density 
of the liquid in atoms per unit volume. Equation 
(4) follows directly from (1), and (3) from the 
assumption that if we take a small imaginary 
sphere of radius R; situated anywhere at random 
in the liquid it will enclose the same number of 
atoms, no matter where it might be in the 
liquid. o may then be found by noting that if in 
Eq. (1) R=R,, then 


y1=3N1/5Ric. (6) 


It is interesting to note that values of o cal- 
culated by this last method are 0.37A for 
T=91.8°K and 0.45A for T=84.3°K (see Table 
Il). These are just the reverse of the values ob- 
tained by the iteration method above. Since the 
R—C proposal depends on the rather precise 
location of Ri and y; from the experimental 
curve, while, on the other hand, by an iteration 


TABLE II. Calculated N’s and o’s (R—C). 














T =91.8°K T =84.3°K 
Ri =3.74A R2=4.6A Ri, =3.75A R2=5.2A 
y1=2.48 4y2= 1.26 91 = 3.36 y2= 1.04 
Ni =5.66 atoms Ni, =9.48 atoms 
N2=3.54 atoms N2=4.07 atoms 
o¢=0.37A o=0.45A 


po=0.0207 atoms/A* po=0.0212 atoms/A® 








0G. S. Rushbrooke and C. A. Coulson, Phys. Rev. 56, 
1216 (1939), : 
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Fic. 2. A comparison of theoretical and experimental 
atomic distribution curves for 4rRp(R) vs. R for liquid 
argon at 91.8°K. Solid line is experimental. Dotted line 
is theoretical (e by R—C method). Circles are theoretical 
values (¢ by iteration). 


process such as has been used, one uses an 
extended portion of the experimental curve, it 
would seem preferable to use the latter method. 
Furthermore, in the cases of liquid sodium and 
potassium,‘ o increased with increasing tem- 
perature which is more reasonable than the 
results obtained by this particular application 
of the R—C method to liquid argon. 

Figures 1 and 2 show theoretical distribution 
functions using the sigmas obtained by both 
methods. Each R —C distribution curve (dotted) 
represents an extended distribution curve and 
hence as indicated in Eq. (2) is a summation, in 
the interval where this is possible, of two 
theoretical distribution equations of the form 
given by Eq. (1); Table II shows the N’s and 
values of ¢ for each curve. These extended R—C 
curves agree reasonably well with the experi- 
mental curves in their general aspects, but in 
both cases the saddle between the experimental 
peaks does not fall to as low a value as does that 
for the theoretical curve. This low value to which 
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the theoretical curve falls in each case is due 
primarily to the fact that the interval in which 
the two theoretical distribution equations add is 
a very small interval, 40 —-(R2—R:). This interval 
in the one case is 0.42A, and in the other is 
0.35A; furthermore, the two functions them- 
selves are quite small in the interval where they 
add. 

The theoretical points represented by circles 
(o by iteration) in Figs. 1 and 2 agree very well 
indeed with the experimental curves and would 
seem to verify the values of o obtained in this 
manner. 


4. THERMODYNAMIC APPLICATION 


Guggenheim" has shown that the free energy 
of a simple liquid is given approximately by the 
equation 


F=—NkT[In J(T)+Inv]+¢, (7) 


where ¢ is the minimum potential energy of the 
molecular assembly and J(T) is the ordinary 
partition function of the molecule. For a mon- 
atomic liquid at moderate temperatures, J(T) 
may be supposed to take its classical value 
[ 2rmkT/h? |! and v is the free volume per atom. 
The other symbols have their usual significance. 

Differentiating Eq. (7), we obtain for the 
entropy of liquid argon 


dF q 
Si= -(—) = wal n J(T)+1n -1o? 
dT/y 3 


1 E. A. Guggenheim, Proc. Roy. Soc. A135, 181 (1932). 
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Using the values og4.3=0.37A and o91.3=0.44A, 
by ordinary linear interpolation’ we have o37.3 
=0.40A for the value of o at the boiling point 
of liquid argon, and 0.07/7.5 =do/dT | rs7.3 for 
an approximate value of do/dT| 787.3. This 
rough interpolation and averaging process is all 
that can be reasonably attempted since the 
temperature range represented by the known 
values of o is quite small. The value of the 
entropy of liquid argon at the boiling point cal- 
culated by this last equation with conversion 
from ergs to calories is 15.8 cal. degree! mole—. 

The approximate entropy of the vapor may be 
calculated from the equation” 


—| V5 


Sy= Ne} in| +In—+-+I1n 2}, (9) 
N 2 


h? 


where V is the volume per mole at the boiling 
point. This calculation gives Sy=32.2 cal. 
degree mole. The entropy of vaporization is 
thus 16.4 cal. degree mole. Since the heat of 
vaporization is 


H=T>).».(Sv—Sz), (10) 
the corresponding heat of vaporization for argon 
at its b.p. is 6.01 kilojoules per mole. The ob- 
served value is 6.3 kilojoules per mole.!* The 
fairly close agreement between the observed and 
calculated values for this heat of vaporization is 
good, considering the approximations used and 
the neglect of communal entropy in the cal- 
culations. 


2 R. H. Fowler, Statistical Mechanics (Cambridge Uni- 
versity Press, New York, 1936), second edition, p. 192, 
Theorem 6.31. 

13 International Critical Tables, Vol. 1. 
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Deviations from Stoichiometric Proportions in Cuprous Iodide* 


RoBert J. MAURER** 
Laboratory for Insulation Research, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received February 23, 1945) 


The absorption of iodine by cuprous iodide has been studied by means of a quartz micro- 
balance. At 132.1°C the concentration of absorbed iodine varies as the square root of the 
pressure of iodine vapor in equilibrium with the sample. The temperature dependence of the 
the reaction at constant iodine vapor pressure indicates a heat of reaction of 0.24 ev per atom 
absorbed. The electrical conductivity of cuprous iodide at 132.1°C varies approximately as the 
4/3 power of the concentration of the absorbed iodine. 
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INTRODUCTION 


FUNDAMENTAL problem in the inves- 

tigation of semi-conducting compounds is 
the dependence of the electrical conductivity 
upon the deviation from stoichiometric propor- 
tions. This deviation is ordinarily so small as to 
prohibit direct measurement by chemical or 
physical methods. 

In 1909, Baedeker! discovered that cuprous 
iodide, at room temperature, will absorb as much 
as 10-* gram of excess iodine per gram of 
cuprous iodide. The resulting electrical con- 
ductivity was electronic in nature but of the 
“defect”’ type, the Hall constant being anomalous 
in sign. Nagel and Wagner® have shown that, if 
cuprous iodide is in equilibrium with iodine 
vapor, at constant temperature, the conductivity 
is a complicated function of the iodine vapor 
pressure. 

In the light of Baedeker’s results it seemed 
feasible to determine the concentration of excess 
























Fic. 1. Quartz microbalance. 





* Publication No. 12 of the Solar Energy Conversion 
Research Program. 

** Now at the Carnegie Institute of Technology, Pitts- 
burgh, Pennsylvania. 

1K. Baedeker, Ann. d. Physik 29, 566 (1909). 

*K. Nagel and C. Wagner, Zeits. f. physik. Chemie 
B25, 71 (1934). 
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iodine in cuprous iodide, as a function of the 
equilibrium iodine pressure and of the tem- 
perature, by direct weighing with a micro- 
balance. Since equilibrium between the cuprous 
iodide and iodine vapor is established in a time 
of the order of minutes, it is necessary to operate 
the microbalance in an atmosphere of iodine 
vapor. | 


EXPERIMENTAL PROCEDURE 


The balance, illustrated in Fig. 1, is an adap- 
tation of the type successfully used for the com- 
parison of gas densities.* 4 It was constructed of 
fused quartz, which was the only material ex- 
posed to the iodine vapor. The balance was 
operated as a null instrument, remote control 
being achieved in the manner described by 
Blewett.® 

The frame, A, was of heavy fused quartz rod. 
The beam, B, was constructed of six pieces of 
half-millimeter quartz rod which were fused 
together, mutually at right angles, with the aid 
of a graphite jig. Two of the pieces formed the 
beam proper and the two short vertical rods were 
used to adjust the center of gravity. The remain- 
ing two quartz rods were cut off short and quartz 
fibers fused to them. These quartz fibers sup- 
ported the beam, being fused to the heavy 
frame. The fibers were slightly slack so as to 
exert negligible restoring torque. The fiducial 
pointer, C, established the zero of the balance. 
The iron cylinder E, which was sealed in a thin- 
walled quartz capillary made an angle of about 


3 A. Stock and G. Ritter, Zeits. f. physik. Chemie A119, 
333 (1926). 

4M. Woodhead and R. Whytlaw-Gray, J. Chem. Soc. 
209, 846 (1933). 

5 J. Blewett, Rev. Sci. Inst. 10, 231 (1939). 
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Fic. 2. Calibration curve of quartz microbalance. 


20 degrees with the beam. This iron cylinder was 
5.0 mm long and 1.0 mm in diameter. 

The horizontal magnetic field of a pair of 
Helmholtz coils was used to create the torque 
upon the iron cylinder which deflected the beam. 
A set of storage batteries supplied the current 
for the coils. 

An important detail in the construction of the 
balance was the method of suspending the speci- 
men, D, and its counterweight, F, from the 
beam. They were hung on quartz fibers which 
were fused to the beam. Only in this manner was 
it possible to maintain the lever arm of the 
balance constant and avoid erratic apparent 
changes in weight of the specimen. 

The counterweight, F, was a piece of platinum 
wire sealed in a quartz capillary. Its density was 
adjusted to within a few percent of the density 
of cuprous iodide so as to avoid buoyancy 
effects. 

The balance was calibrated by placing known 
weights of aluminum wire on the specimen hook 
and determining the current necessary to bring 
the beam back to its zero. The aluminum wire 
was weighed upon an Ainsworth microbalance. 
The calibration was performed by first placing 
upon the specimen hook an unknown mass suf- 
ficient to bring the beam pointer slightly below 
the fiducial pointer. The current to balance the 
beam was %. The known mass, m, was then 
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added to the specimen hook and the current for 
balance was 7. Figure 2 shows a calibration 
curve. In agreement with the theory of the 
balance, 7?—79? was proportional to m. The zero 
of the balance was observed over a period of 
several weeks and remained constant within a 
fraction of a microgram. 

The reliability and precision of the instrument 
were tested by mounting it in a Pyrex tube with 
a molybdenum weight on the specimen hook. At 
room temperature, the apparent change in 
weight of the molybdenum was determined as a 
function of the pressure of dry nitrogen in the 
tube. The data are shown in Fig. 3. 

Figure 4 shows the experimental arrangement 
for determining the dependence of the weight of 
a cuprous-iodide specimen upon the equilibrium 
pressure of iodine vapor. The balance was in the 
Pyrex tube, JT. The Helmholtz coils H, sur- 
rounded the Pyrex tube. The entire assembly 
was in the electrically heated oven, O. Two 
small, double-walled, mica windows permitted 
observation of the beam and fiducial pointer. 
A simple lens giving a magnification of 5-fold 
was used for this purpose. The Pyrex.tube, E, 
which connected the balance tube with the iodine 
reservoir, A, was wrapped with a heating element 
and kept at a temperature intermediate of that 
of the balance and iodine. 

Before assembly of the apparatus a sample of 
0.5-mil copper foil was suspended from the 
specimen hook. The weight of foil chosen was 
such as to yield a weight of cuprous iodide ap- 
proximately 10° gram in excess of that necessary 
to balance the counterweight. The capsule, D, 
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Fic. 3. Test of quartz microbalance 
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contained iodine which had been presublimed in 
a stream of dry nitrogen. After evacuating the 
tube through B, the capsule D was opened by 
dropping the glass-enclosed iron hammer S, on 
the thin-walled glass tip, C. The iodine was then 
distilled into the lower part of the tube E. 
Finally, after pumping the tube to a pressure of 
10-* mm Hg, it was sealed off at A. 

The copper foil was slowly transformed to 
cuprous iodide by holding the balance oven at 
150°C and the iodine oven at 80°C. 

The oven temperatures were held constant by 
mercury thermoregulators and relays and the 
vapor pressure of the iodine in the tube was 
determined from the temperature of the iodine 
oven by use of the vapor pressure data of the 
International Critical Tables. 

While taking data, the zero of the balance was 
checked periodically by immersing the iodine 
in liquid air and thus reducing the vapor pressure 
of iodine in the tube to a negligible value. 

Figure 5 shows the dependence of the density 
of absorbed iodine atoms upon the pressure of 
the iodine vapor in equilibrium with the cuprous 
iodide. The temperature of the cuprous iodide 
was 132.1°C. The density of the absorbed iodine 
was calculated from the initial weight of the 
pure cuprous iodide and its measured change in 
weight. Since the initial weight of the copper foil 
was 29.414 milligrams, the initial weight of the 
cuprous iodide was taken as 88.14 milligrams. At 
all but the lowest iodine vapor pressures, the 
density of iodine atoms is proportional to the 
































Fic. 4. Apparatus for investigation of absorption of iodine 
by cuprous iodide. 
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square root of the iodine vapor pressure. 
d= dopte}. (1) 


Figure 6 shows the dependence of the density 
of absorbed iodine atoms upon the temperature 
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Fic. 5. The density of absorbed iodine atoms as a func- 
tion of the equilibrium pressure of iodine vapor at constant 
temperature, ¢=132.1 degrees C. 


of the cuprous iodide for a constant equilibrium 
vapor pressure of iodine of 1.55 mm Hg. Over 
the narrow temperature interval examined, the 
density of absorbed iodine atoms is an exponen- 
tial function of the reciprocal of the absolute 
temperature. The slope of the logio d versus 1/T 
curve shown in Fig. 6 gives the heat of reaction 
as 0.24 electron volts per atom absorbed. 


ELECTRICAL CONDUCTIVITY 


The electrical conductivity of cuprous iodide 
was also measured as a function of temperature 
and of the equilibrium pressure of iodine vapor. 
Nagel and Wagner have shown that the in- 
creased conductivity, which results from the 
absorption of iodine, is electronic in nature. The 
electrolytic conductivity (due to the Cut) is also 
increased but is negligible. The experimental 
arrangement, which was essentially that of 
Nagel and Wagner,” is illustrated in Fig. 7. 
Platinum-foil electrodes were spot welded to 
two strips of copper foil, whose dimensions were 
2.0 cm X0.5 cm X 1.2K 10-* cm. The copper foils 
were suspended on platinum springs, spot 
welded to the electrodes, inside a sealed-off, 
evacuated Pyrex tube. Iodine, which had been 
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F Fic. 6. The density of absorbed iodine atoms as a 
function of the temperature at constant iodine pressure. 
pl2=1.55 mm Hg. Heat of reaction=0.24 ev per atom 
absorbed. 


presublimed in a stream of dry nitrogen, was 
distilled into the tube prior to its final evacuation 
and sealing off. That part of the tube containing 
the copper foils was surrounded by an electrically 
heated oven, O;. The tube, E, which connected 
the specimen chamber, S, with the iodine reser- 
voir, D, was wrapped with an electrical heating 
element and kept at a temperature intermediate 
of that of the ovens, O;, and O2. The oven tem- 
peratures were controlled by thermoregulators 
and measured with thermocouples and mercury 
thermometers. The resistance of the cuprous 
iodide was measured with a Leeds and Northrup 
Wheatstone bridge. The independence of the 
resistance of the cuprous iodide of the magnitude 
and sign of the d.c. voltage across the sample 
was checked. The resistance as measured with 
60-cycle alternating current was the same as the 
d.c. resistance, within the precision of measure- 
ment. 

The data are collected in Fig. 8. The agree- 
ment between these results and those of Nagel 
and Wagner is good except at the lowest iodine 
pressures where Nagel and Wagner’s conduc- 
tivities are as much as 50 percent lower than 
those found in these experiments. 

By interpolation, the conductivity of cuprous 
iodide as a function of iodine pressure for a 


MAURER 

cuprous-iodide temperature of 132.1°C, can be 
obtained from the data of Fig. 8. In conjunction 
with the results of Fig. 5, the conductivity as a 
function of the density of absorbed iodine atoms, 
can be plotted. This is done in Fig. 9. Over the 
range of excess iodine concentration investigated, 
the conductivity varies approximately as the 4/3 
power of the density of the absorbed iodine 
atoms. 


o =o 04/8, (2) 


DISCUSSION 


Nagel and Wagner have proposed a mechanism 
for the absorption of iodine by cuprous iodide, 
similar to that which has been used, to success- 
fully describe the behavior of cuprous oxide. 
They assume that two cuprous ions and two 
electrons wander to the surface of the solid for 
each iodine molecule absorbed. The result is an 
extension of the cuprous-iodide lattice, the 
formation of two vacancies in the cuprous ion 
lattice, and a defect of two electrons. This 
mechanism is plausible because the slight elec- 
trolytic conductivity of cuprous iodide is entirely 
due to the mobility of the Cut ion. The reaction 
can be described by the chemical equation 


2 Cut+1.(g)=—2 Cut+2I-+2V+2h, (3) 


where the symbol, V, denotes a cuprous ion 
vacancy and the symbol, hf, a defect of an 
electron or a ‘‘hole.”” The concentration of 
‘‘holes” will then be proportional to the fourth 
root of the iodine vapor pressure, if the ‘‘holes’’ 


























Fic. 7. Apparatus for determination of the electrical 
conductivity of cuprous iodide as a function of temperature 
and equilibrium vapor pressure of iodine. 
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are assumed free. 
(h) =n = nop I". (4) 
Nagel and Wagner assume that the conductivity 
is proportional to the concentration of ‘‘holes,”’ or 
o=dopl 2°. (S) 


This result is not in agreement with the experi- 
mental data of Fig. 8, although there is some 
indication that when the equilibrium vapor 
pressure is very small, Eq. (5) may be valid. 
The conductivity of a semi-conductor is 
usually discussed in terms of Lorentz’s equation: 


o = (4/3)[e*nl/(2rmkT)*]. (6) 


The conductivity is proportional, at constant 
temperature, to the product of the density of 
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Fic. 8. The electrical conductivity of cuprous iodide as 
a function of temperature and equilibrium iodine pres- 
sure. 
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Fic. 9, The electrical conductivity of cuprous iodide as a 
function of the density of absorbed iodine atoms. 
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Fic. 10. The mean free path of the conducting “‘holes’’ 
in cuprous iodide as a function of the density of “‘holes.” 
After Steinberg. 


charge carriers, n, and their mean free path, /. 
In Fig. 10, is shown the dependence of the mean 
free path, /, of the ‘‘holes” in cuprous iodide, 
upon the density of ‘“‘holes.’’ These data are 
calculated from the Hall effect and conductivity 
data of Steinberg.® If we assume that the density 
of ‘‘holes” is of the order of magnitude of the 
density of absorbed iodine atoms, the mean free 
path is a strong function of the density of 
‘tholes’”’ for the range of densities involved in 
these experiments. In fact, the mean free path 
can be well respresented within this range, by 
the equation 

L=lon—®-3, (7) 


If we assume that this relationship, obtained 
from Steinberg’s measurements at room tem- 
perature, is valid at the higher temperature of 
the present experiments, we can combine (7) and 
(6) and conclude that 


o = (4/3)[e%lon®-7/(2amkT)*], (8) 


combining Eqs. (2) and (8), we have the result 
that, within the error of experiment, the density 
of conducting “‘holes’’ varies as the square of 
the density of absorbed iodine atoms, 


nx d*; (9) 


and is proportional to the equilibrium vapor 
pressure of iodine. 
n « plo. (10) 


6 K. Steinberg, Ann. d. Physik 35, 1009 (1911). 
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CONCLUSION 


The effect of the absorption of iodine upon the 
electrical conductivity of cuprous iodide is not 
adequately explained by the mechanism of 
Nagel and Wagner. We plan to extend the ab- 
sorption measurements to other temperatures 


C. HOWARD SHOMATE 


and to determine the Hall coefficient of cuprous 
iodide as a function of temperature and concen- 
tration of absorbed iodine. 

The author is deeply indebted to Professor A. 
von Hippel, whose advice and assistance made 
possible the performance of the experiments 
described in this paper. 
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The specific heats of a- and y-manganese were measured in the temperature range 52° to 
298°K. The specific heat curve of a-manganese exhibited a small ‘‘hump,” extending over the 
temperature range 63°-103°K, with the peak of the “hump” occurring about 95°K. No 
anomaly appeared in the specific heat curve of y-manganese. The computed heat contents and 
entropies per gram-atom are as follows: a-manganese, H°29s.16—H°9=1193+4 cal., and 
S°298.16= 7.59+0.04 E.U.; y-manganese, H°20s.16s—H°o=1221+44 cal., and S°293 16=7.72 


+0.04 E.U. 


HIGH temperature allotropic modification 

of a metal seldom can be obtained in a suf- 
ficiently pure and stable form to determine its low 
temperature specific heat. Only one instance has 
been reported previously in the literature; viz., 
gray tin.? y-manganese which is stable at low 
temperatures has been prepared by a controlled 
electrolysis at the Salt Lake City Station of the 
Bureau of Mines. This paper presents low tem- 
perature specific heat data for this material and 
for a-manganese. 


MATERIALS 


The y-manganese used in these measurements 
was prepared by David Schlain‘ by electrolysis 
of manganous sulfate in the absence of sulfites. 
In the presence of sulfites the electrolysis de- 
posits a-manganese. The resistivity of each of 
forty-three strips of y-manganese used in the 
measurements was determined by E. V. Potter.5 
The mean value was 45.1 microhms-cm, with a 

1 Published by permission of the Director, Bureau of 
Mines, U. S. Department of the Interior. 

2 Chemist, Pacific Experiment Station, Bureau of Mines. 

3F. Lange, Zeits. f. physik. Chemie 110, 343 (1924). 


‘ Metallurgist, Salt Lake City Station, Bureau of Mines. 
5 Physicist, Salt Lake City Station, Bureau of Mines. 


standard deviation of 0.4 microhms-cm. Previous 
measurements of other samples of y-manganese 
showed resistivities of 42.8+0.3, 45.6+0.5, and 
44.7+0.2 microhms-cm. The resistivity of a-man- 
ganese is 185 microhms-cm. This is rather con- 
clusive evidence of the absence of a-manganese 
in the y-manganese sample; certainly it is not 
present to the extent of a few percent. The 
sample of y-manganese was refrigerated im- 
mediately after preparation and, except for a 
few hours at room temperature, was kept re- 
frigerated until the specific heat measurements 
were commenced. Indications are that the rate 
of conversion of y- to a-manganese is very slow, 
even at room temperature. A 374.20-g sample 
of y-manganese was used in the measurements. 

After completion of the measurements on 
y-manganese, the sample was removed from the 
calorimeter and conversion to a-manganese 
effected by heating gradually to 550°C in a high 
vacuum (10-‘ mm Hg). Before conversion, the 
y-manganese was ductile, but after conversion 
to a-manganese it became brittle. Specific heat 
measurements were then made on a 368.30-g 
sample of this a-manganese. 
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SPECIFIC HEATS OF a-MANGANESE AND y-MANGANESE 




































































ous as 0.5 percent or as high as 2.0 percent, the cor- 
en- nn rected specific heat results would still be essen- 
Oo 6 ee oaae tially the same as those obtained assuming’ 1.0 
A, W ALB percent impurities. The corrected specific heat 
ade « Pa values of y-manganese are also listed in Table I. 
nts y 
lu / SPECIFIC HEATS 
= * The method and apparatus used in the low 
2 Fd temperature specific heat measurements have 
* been described previously.* * The experimental 
45 a f results in Table I are expressed in defined calories 
z 2 (1 calorie=4.1833 international joules).!° The 
O f results are represented graphically in Fig. 1. 
a ] Values of the specific heats at 298.16°K, ob- 
Vv tained by extrapolation of smooth curves 
through the experimental .points, also are in- 
0 100 S00 Sto cluded in Table I. The atomic weights are in 
™%* accordance with the 1941 International Atomic 
: Weights. All weights were corrected to vacuum. 
Fic. 1. Specific heats of y-manganese (A), and : 
conan Oy: TABLE I. Specific heats of manganese 
(at. wt. = 54.93 g). 
The dissolution of the sample of y-manganese 
in 1N H2SO« left a residue of 0.4 percent. In rm g Manganese | som . Cae 
addition, it was possible to detect a slight odor 
us of acetylenic hydrocarbons during the dissolu- -s ey 74 = 
se tion. It was felt that the best. way to make cor- 62.1 1.806 61.7 1.829 
id rections for these impurities was to make specific fe ia Te po 
n- heat measurements of an available sample of very 82.2 2.928 76.1 2.555 
oll pure, electrolytic a-manganese prepared at the oe oe oy po 
se Bureau of Mines Station at Boulder City, 93.3 3.563 94.7 3.340 
ot Nevada.*? The purity of this sample is pur- 94 ane TE = 
ne ported to be at least 99.9 percent. It was degassed 101.2 3.538 124.3 4.293 
ae at 750°C in a high vacuum (10-‘ mm Hg): The tg oon ry oa 
a sample of this material used in the measure- 129.3 4.336 155.5 4.979 
c- ments weighed 484.61 g, and the specific heats oy os a7 ase 
ts are tabulated in Table I. The specific heat 190.4 5.399 185.6 5.458 
fe results for y-manganese were corrected for im- oy re ay eae 
, purities by assuming 1 percent impurities and 241.2 5 886 216.5 5.851 
le correcting the results by a comparison of the os re Soy oe 
S. specific heats of a-manganese prepared by con- 271.0 6.119 246.1 64.131 
n 7 version of y-manganese and the Boulder City ap erry ee a. 
e . a-manganese. The corrections raise the specific 301.3 6.301 276.4 6.414 
e heat values of the uncorrected y-manganese (298.16) (6.29) et zo 
h about 0.2 percent on the average. If it had been (298.16) (6.59) 
e assumed instead that the impurities were as low 
a ®R. S. Dean, Metals Technology 11, No. 4, Technical 8K. K. Kelley, J. Am. Chem. Soc. 63, 1137 (1941). 
it > Paper No. 1721 (1944), 17 pp. *C. H. Shomate and K. K. Kelley, J. Am. Chem. Soc. 66, 
g "J. H. Jacobs, J. W. Hunter, W. H. Yarroll, P. E. 1490 (1944). 4a 
Churchward, and R. G. Knickerbocker, Metals Technology 10 EF. F. Mueller and F. D. Rossini, Am. J. Phys. 12, 1 


11, No. 5, Technical Paper No. 1717 (1944), 21 pp. (1944). 
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TABLE I]. Heat contents and entropies at 298.16°K. 








a-Manganese 


E.U./g-atom 


H° —H,°, 
cal./g-atom 


—" y-Manganese 
—H°, 


cal./g-atom E.U./g-atom 





18.3 
1175. 
1193.44 


0°-52°K (extrap.) 
52°-298.16°K (measured) 
0°-298.16°K 


0.49 
7.23 
7.72+0.04 


18.9 
1202. 
1221.44 


0.48 
7.11 
7.59+0.04 











The specific heat curve of y-manganese showed 
no anomalies in the temperature range studied. 
However, a small ‘‘hump”’ appears in the spe- 
cific heat curve of a-manganese, extending over 
the temperature range 63°-103°K, with the peak 
of the ‘‘hump”’ occurring about 95°K. 

Kelley" reported low temperature specific heat 
data for a-manganese, which had been prepared 
by electrolyzing a solution containing manganous 
sulfate and ammonium sulfate into a mercury 
cathode and distilling the mercury from the 
resulting amalgam. The sample was a grayish 
powder consisting of clumps of particles of 4- to 
5-microns average size. The low apparent density 
resulted in a sample size of only 112.53 g, which 
is less than one-fourth the size of the sample in 
the present measurements. In addition, oxide 
impurities corresponding to about 1 percent 
MnO were present. With these considerations the 
agreement of Kelley’s work with the present 
measurements may be considered good. In 
general his results were about 0.3 percent higher. 
A few of Kelley’s representative values have been 
plotted (black triangles) in Fig. 1. 

Kelley did not report a “hump” in his specific 
heat curve of a-manganese. However, recalcula- 
tion of his original data showéd his specific heat 
value at 88.2°K to be in error. The corrected 
value is 3.251 cal./g-atom, or about 4 percent 
higher than the value of 3.130 cal./g-atom which 
he reported. Thus, re-examination of his specific 
heat curve indicated a “thump” of magnitude 
comparable to that in the present measurements. 


HEAT CONTENTS AND ENTROPIES 


The evaluation of the heat contents at 298.16°K 
is made by graphical integration from a plot of 
C, against 7, and entropies from a plot of C, 
against log T. This necessitates extrapolation of 


1K. K. Kelley, J. Am. Chem. Soc. 61, 203 (1939). 


the specific heat curves from the temperatures 
of the lowest measurements down to the absolute 
zero of temperature. The specific heat curves of 
a-manganese and -manganese were each fitted 
at the lower temperatures of the measurements 
by a single Debye function. 

The values of @ in the Debye function, D(6/T), 
were calculated for the a-manganese specific heat 
determinations at 54.2°, 58.2°, 62.1°, and 104.3°K 
to be 347.4, 348.0, 347.7, and 347.4, respectively. 
The values for intermediate determinations could 
not be used, since these specific heat values are 
affected by the ‘‘hump.’’ The mean value is used 
to extrapolate the specific heat curve of a-man- 
ganese by the function, D(347.6/T), yielding the 
extrapolated portion of the heat content, 
F°52—H°y=18.3 cal./g-atom, and of the en- 
tropy, S°s2=0.48 E.U./g-atom. 

Evaluation of @ from the nine lowest tem- 
perature specific heat determinations of y-man- 
ganese gave a mean value of @= 343.5+1.0. The 
function, D(343.5/T), gives the extrapolated 
values H°5.—H°)=18.9 cal./g-atom, and S°52 
=0.49 E.U./g-atom. Above 100°K the specific 
heat values of both a-manganese and y-man- 
ganese begin to deviate from the Debye functions 
used in the extrapolations, the deviations becom- 
ing quite marked at higher temperatures. 

The heat content and entropy associated with 
the “hump” in a-manganese was determined by 
summing the energies of eight successive specific 
heat determinations over the temperature range 
80° to 106°K in the manner described by 
Shomate.” The excess energy in the ‘‘hump,” 
above the specific heat curve represented by the 
function D(347.6/T), is 4.72 cal./g-atom, and 
the excess entropy 0.054 E.U./g-atom. 

Table Il summarizes the heat content and 
entropy calculations for a- and y-manganese. 


2 C. H. Shomate, Ind. Eng. Chem. 36, 910 (1944). 
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Heat contents above 298.16°K of two samples of electrolytic manganese were determined 
from room temperature to about 1440°K. From these data three transitions occurring at 
1000+2°, 1374+5°, and 1410+5°K with heat effects of 535, 545, and 430 calories per gram- 
atom, respectively, were observed. Smooth curve values at 100° intervals and corresponding 
entropies are given as well as heat content and specific heat equations for the four modi- 
fications of manganese. A specific heat equation for y-manganese was derived for the tem- 
perature range 298-1374°K from low temperature and high temperature thermal data; and 
corresponding heat contents and entropies above 298.16° for y-manganese were calculated. 



















producing erratic values. Heat content measure- 


HE heat content of manganese at high 
ments, recently completed at this station, of two 


temperatures and its transition points have 





been studied by numerous investigators,* but 
their data are inconsistent. This disagreement 
: would appear to be due largely to lack of sample 
d purity, resulting in sluggish transitions and thus 







TABLE I. Heat content of manganese (at. wt. = 54.93). 











T —H 293.16 Hr —H293.16 


Hr —H 293.16 H 
T, °K (cal./g-atom) 
































T,°K (cal./g-atom) T,°K  (cal./g-atom) 
Sample A 

500.9 1406 993.5 5385 1354.2 9200 
555.8 1822 999.2 5880 1355.6 9240 
634.7 2378 1002.0 5980 1371.3 9595 
648.8 2483 1025.9 6230 1373.6 9880 
658.8 2563 1056.8 6490 1388.7 10,120 
758.5 3374 1111.6 7010 1394.0 10,155 
870.9 4328 1153.3. 7375 1405.0 10,250 
927.9 4809 1220.4 7975 1413.7 10,800 
940.0 4917 1281.5 8525 1416.1 10,840 
956.2 5060 1301.5 8725 1427.3 10,960 

4 962.9 5125 1321.2 8945 1436.3 11,060 
986.0 5300 1336.4 9075 

Sample B 

: 366.9 469 1000.8 5605 1386.2 10,060 

f 520.1 1558 1003.0 5990 1386.3 10,000 
650.5 2515 1074.6 6620 1410.3 10,645 
802.4 3683 1291.9 8585 1410.6 10,510 
822.8 3867 1355.2 9260 1411.6 10,400 
970.0 5155 1369.8 9490 1414.1 10,320 
990.3 5330 1376.5 9840 1418.7 10,505 
996.3 5485 1382.1 9950 














1 Published by permission of the Director, Bureau of 
Mines, U. S. Department of the Interior. 

* Chemist, Pacific Experiment Station, Bureau of Mines. 

3 (a) J. C. Southard and C. H. Shomate, J. Am. Chem. 
Soc. 64, 1770 (1942); (b) A. Sieverts and H. Moritz, Zeits. 
f. physik. Chemie A180, 249 (1937); (c) H. Yoshisaki, 
Sci. Rep. Tohoku Imp. Univ. 26, 182 (1937); (d) L. V. 
Gayler, J. Iron Steel Inst. 115, No. 1, 393 (1927); (e) S. 
Umino, Sci. Rep. Tohoku Imp. Univ. 16, 775 (1927); 
(f) F, Wiist, A. Meuthen, and R. Durrer, Forsch. Arb. 
Ver. deut. Ing. No. 204 (1918); (g) R. Laemmel, Ann. d. 
Physik 16, 551 (1905); (hk) N. Stucker, Sitz. kénig. Akad. 
Wiss. Wien 114, 657 (1905). 
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samples of pure electrolytic manganese gave 
results in good mutual agreement, from which 
the temperature and thermal magnitude of three 
transitions were determined. 


METHOD AND MATERIAL 


The heat content measurements were made by 
the “‘drop’’ method in an apparatus previously 
described. Frequent calibrations of the furnace 
thermocouple were made at the gold point in the 
manner outlined by Southard. The heat capacity 
of the copper calorimeter was determined elec- 
trically, the results being converted to the ther- 
mochemical calorie by the relation,’ 1 cal. 
= 4.1833 int. joules (VBS). During the measure- 
ments the samples were enclosed in fused silica 
capsules filled with helium before sealing. The 
heat content of the capsule material was deter- 
mined previous to the manganese measurements 
and closely checked values already reported* ® 
for silica. Because of rapid devitrification of the 
silica above 900°C about 10 capsules were used 
during the course of the measurements and 
filled with an unused portion of manganese each 
time. The heat content of the silica contributed 
about 35 to 40 percent of the total heat evolved 
in the calorimeter. 

The two electrolytic manganese samples were 
prepared at the Boulder City Station of the 


4J. C. Southard, J. Am. Chem. Soc. 63, 3142 (1941). 

5 E, F. Mueller and F. D. Rossini, Am. J. Phys. 12, 1-7 
(1944). 

6W. P. White, Am. J. Sci. 47, 1 (1919); J. Am. Chem. 
Soc. 55, 1050 (1933). 
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Fic. 1. High temperature heat contents of manganese 
above 298.16°K. 


Bureau of Mines and were at least 99.9 percent 
pure after degassing free of hydrogen. Sample A 
was degassed by heating to 850°C in a vacuum 
of 10-5 mm and allowed to cool slowly. Sample B 
was treated similarly at 900°C. Although the 
impurities present in each sample are very small, 
sample A was considered the purer, as it was 
free of silicon while sample B contained 0.03 
percent. 


RESULTS 


Experimental heat contents above 298.16°K 
are given in Table I for both samples A and B. 
Except for several low points in the 6-range, the 
results of the two samples are in good agreement. 
From the data, the existence of three transition 
points in the range 298°-1436°K is readily dis- 
cernible. These transitions, as well as the experi- 
mental points, are shown in Fig. 1. The listing 
arrangement in Table I does not necessarily 
conform with the experimental order; in fact, 
many of the determinations at lower tempera- 
tures followed those at higher temperatures. 

It was assumed that the manganese came to 
the same final state (alpha) in the calorimeter, 
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regardless of its temperature before cooling. For 
all determinations in Table I within 50° of the 
a—8 transition point, the furnace was brought 
to thermal equilibrium; the capsule was then 
admitted and, upon reaching furnace tempera- 
ture, was dropped into the calorimeter. Thus the 
capsule always was heating before dropping. 
This procedure was followed upon finding that 
6-manganese could be readily supercooled. Such 
sluggishness was not noticed in the transitions at 
higher temperatures. 

Series of measurements, made to observe the 
extent to which 6-manganese could be super- 
cooled, are arranged in Table II in the order of 
determination. These involved heating the 
sample to at least 1010°K within the furnace, 
cooling the furnace as rapidly as possible to the 
run temperature, bringing it to equilibrium, and 
then proceeding to drop the sample as usual. The 
cooling process required about 20 to 30 minutes. 
It was found that 6-manganese could be super- 
cooled to 960°+10°K. Under faster cooling con- 
ditions a reduction of this value might be 
achieved. 

The transition points and associated heat 
effects obtained from the data given in Table I 
are tabulated in Table III. 

Smooth curve values at 100° intervals and cor- 
responding graphically computed entropies are 
listed in Table IV. Heat content equations repre- 
senting the experimental data have been derived 
and are given below. They are based on the fol- 


TABLE II. Heat content of supercooled 6-manganese 
(at. wt. = 54.93). 


























Ar —H 293.16 Ar —H 2098.16 
7. (cal./g-atom) \ i 4 (cal./g-atom) 
Sample A Sample B 
972.9 5770 997.2 5880 
961.8 5640 994.2 5915 
950.6 5030 987.4 5860 
979.2 5770 
967.6 5295 
963.0 5520 
TABLE III. 
Heat of trans. 
Transition 7 (cal./g-atom) 
a—>B 1000+2 535 
Boy 1374+5 545 
y> 5 1410+5 430 
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lowing values: a-manganese, Cy 298.16=6.29,7 
Hr—H298.16=2895 and 5450 cal. at 700° and 
1000°K ; B-manganese, Hy —Ho95,1¢ = 5985, 7795, 
and 9395 cal. at 1000°, 1200°, and 1374°K; 
y-manganese, Hr=Ho3.16=9940 and 10,330 
cal. at 1374° and 1410°K; 6-manganese, Hr 
— Ho9s.16= 10,760 and 11,210 cal. at 1410° and 
1450°K, respectively. 

The appropriate temperature range and the 
mean percentage deviation of the equation from 
the experimental data are given after each 
equation. 


Mn(a) : Hr — Ho8.15=5.707 +0.00169T? 
+ (37,500/T) — 1975(298°-1000°K ; 1%), 
Mn(8) : A 7—H_98.16 = 8.33T +0.00033T? 
— 2675(1000°-1374°K ; 0.2%), 
Mn(y7) ° Hr—Hags.16= 10.70T 
—4760(1374-1410°K ; 0.1%), 
Mn(6) ; Hr—H23.16= 11.30T 
—5170(1410°-1450°K ; 0.0%). 


Differentiating the above heat content equa- 
tions gives the corresponding’ specific heat rela- 
tionships. 


Mn(a): C,=5.70+0.003387 — (37,500/T?), 
Mn(8): Cp=8.33-+0.00066T, 

Mn(y7) : C,=10.70, 

Mn(6): C,=11.30. 


DISCUSSION 


Heat content data for manganese comparable 
to this work were reported by Southard** and by 


TABLE IV. Heat content and entropy. 











Hr —H293.16 Sr —S298.16 
7. =< (cal./g-atom) (cal./deg./g-atom) 
400 690 1.99 
500 1385 3.54 
600 2120 4.88 
700 2895 6.07 
800 3715 7.16 
900 4570 8.17 
1000 5450 (a) 9.10 
1000 5985 (6) 9.63 
1100 6890 _ 10.50 
1200 7795 11.28 
1300 8715 12.02 
1374 9395 (8) 12.53 
1374 9940 (vy) 12.93 
1400 10,220 13.13 
1410 10,330 (y) 13.21 
1410 10,760 (6) 13.51 
1450 11,210 13.83 











7C. H. Shomate, J. Chem. Phys. 13, 326 (1945). 


TABLE V. ‘Heat content and entropy of y-manganese. 











Ar —H23.16 Sr —S298.16 
rT, “K (cal./g-atom) (cal./deg./g-atom) 
400 705 2.02 
500 1445 3.68 
600 2230 5.10 
700 3050 6.37 
800 3915 7.52 
900 4815 8.58 
1000 5750 9.57 
1100 6720 10.50 
1200 7730 11.38 
1300 8780 12.21 
1374 9570 12.80 








Wiist and co-workers.® Southard’s results are 
somewhat lower except around 650°K where they 
are in good agreement. From 800°K to the a8 
transition the difference is about 2 percent, thus 
reducing the heat of transition in the present 
instance to 535 calories. In the 8-range there is 
fair agreement if the 65-calorie correction, ad- 
vocated by Southard, is added to his values. 
The values of Wiist are lower than the present 
work by only 1 to 2 percent below 1000°K but 
are in complete disagreement above that tem- 
perature, since no perceptible a—§ transition 
was found. 

The present data for manganese are con- 
sidered more reliable and complete than any 
previously reported. Besides having samples of 
exceptional purity, several refinements in the 
apparatus were made since Southard’s work, and 
considerably more effort was expended to isolate 
and investigate the transitions, especially those 
occurring at higher temperatures. That a fourth 
modification (delta) of manganese does exist is 
shown conclusively by the present data. 

The herein reported transition temperatures 
have been checked within a few degrees at the 
Salt Lake City Station of the Bureau of Mines by 
thermal-analysis and hydrogen-solubility meas- 
urements. The a— transition also was observed 
to be sluggish, its temperature depending on 
whether the sample was heating or cooling at the 
time of measurement. Hence, some doubt exists 
as to the true equilibrium temperature. However, 
for purposes of making thermodynamic caicula- 
tions, 1000°K was selected ; admittedly this may 
produce some error, but it is negligible except 
near or at the transition temperature. 











332 JOHN E. 

Gamma-manganese is found often as a separate 
phase in manganese alloys, but is thermody- 
namically unstable in the pure state below 
1374°K. However, low temperature specific heats 
have been determined,’ as at temperatures less 
than 273°K its rate of transformation is very 
low. Transformation occurs almost instanta- 
neously around 398°K; hence no heat contents 
can be measured between room temperature and 
1374°K. Knowing the total virtual entropy of 
-manganese at 1374°K to be the same regardless 
of its origin, a specific heat equation for y-man- 
ganese was derived, using the data: 


(y-Mn), Ce 298.16 = 6.59, (dC,/dT) 298.16 = 0.0069, 


WALTER 


and 


Soo98.16= 7.72; (a-Mn), S298.16 = 7.59; 


and 
Sisza(y-Mn) — S29s.16(a-Mn) = 12.93, 
C,=6.03+0.00356T — (44,300/T?). 


Integrating gives 


Hr —H2s,16=6.03T +0.001787? 
+ (44,300/T) — 2105, 


where Hr—H293.16 represents the heat liberated 
in calories per gram-atom in cooling y-manganese 
from T°K to produce y-manganese at 298.16°K. 
Heat content values, calculated from this equa- 
tion, and entropies, obtained by integration, are 
given in Table V at 100° intervals to 1374°K. 
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The process of chromatographic adsorption is studied from the kinetic viewpoint, assuming 
a bimolecular reaction between adsorbent and solute. Boundaries characteristic of equilibrium 
adsorption are obtained only if the time of passage of the solution through the column is 
greater than the half-life of the reaction by a factor of at least one hundred. 


INTRODUCTION 


N a recent article! the author presented the 
results to be expected if one or more solutes 
are adsorbed by a chromatographic process, 
assuming the adsorption isotherm to be the 
equilibrium equation for a*bimolecular chemical 
reaction. The distribution of adsorbed solute 
along the adsorbing column was found to be a 
definite function of the equilibrium constants of 
the reactions. In this article the same process is 
investigated from the kinetic viewpoint in order 
to see how the previous results are modified if 
equilibrium is not attained. 
We assume the adsorption to be represented 
by a bimolecular reaction of the type 


HZ +MX=MZ+ 
(P-Q) (C) (Q) 


1 J. Walter, J. Chem. Phys. 13, 229 (1945). 


mA 4 
(C°—C) 


(1) 


where (C) is the concentration of the solute in 
solution, (Q) is the concentration of adsorbed 
solute, (P) is the capacity per unit length of the 
tube, and (C°) is the original concentration of the 
solute. If / is the distance down the tube and v 
is the volume of solution which has passed the 
point /, we have 


80/dv+aC/al=0, (2) 


as the equation of continuity. If @ is the time 
during which the solution has been in contact 
with the adsorbing material at a given position 
along the tube, we have 


0Q/d0=k;C(P—Q) —k,Q(C°—C), (3) 


where ky and ky are the rate constants for the 
forward and backward reactions in (1). If R is 
the rate at which the solution is passed through 
the tube, we have 


80/80 = (dv/d0)(8Q/dv) = RdQ/dv, (4) 
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and Eq. (3) becomes 

RdQ/dv=k;C(P—Q) —kiQ(C’—C). (5) 
Placing R=0 gives the equilibrium relation 
between C and Q, and the problem reduces to 


that previously treated. Eliminating C from 
Eqs. (2) and (5) leads to 


[k;P+(ko—ks)Q ?(0Q/dv) +kskyC°P(9Q/dl) 
rag 


+R Ce,P+(bs—k,) 0} 
dvol 


00 ~0 
ate = =0, (6) 
dv al 


as the differential equation governing the dis- 
tribution of Q along the adsorbing column. The 
boundary conditions are Q=0 for v=0 and 
C=C at 1=0. This last condition implies 


RdQ/dv=k;C(P—Q); 1=0, vZ0, (7) 
Q=P{1-—exp [—k;C%/R]}; 1=0. (8) 


As a subsidiary condition which can be used as 
a check, we have . 


ff ou- C%, (9) 


expressing the conservation of solute. The sub- 
stitutions 


q=Q/P, 


reduce the equations to 


v=Rt/Ck;, 1=Rx/Pk;, (10) 


fe) 
[ky +(ko—kys)q 





q dq 
+koks— 
ot Ox 


9 


“qd 
+hyLky+ (ks —ks)g }— 
oxdt 


dq 9q 
—k;(ky—k;)— —=0; (11) 
Ot Ox 
q=1-—e'; x=0; (12) 
ff aae=e. (13) 
0 


Equation (11) can be solved exactly for k,=0 
and for k,=k,;, although the latter solution 
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cannot be given in closed form. The author has 
not succeeded in obtaining an exact solution for 
the general problem, but for the interesting case 
k,>k; an approximate solution, valid for large 
values of ¢, is presented. 


k,=0 


For k,=0, the term in (11) involving dg/dx 
drops out of the equation, which can then be 
rearranged to give 


re) 1 oq 
a-o—|a+——|=0, (ag 
ot 1—q 0x 


so that either g=1 or 
q+[1/(1—g) ](0q/dx) =f(x), (15) 


where f(x) is an arbitrary function of x. The 
boundary conditions can be satisfied by placing 
f(x) =0. We then have 


q/(1—g) =g(te*, (16) 


where g(t) is an arbitrary function of ¢. From 
(12) we find 


q/(1—g) = (e'—1)e*. (17) 


This solution satisfies Eq. (13), and thus justifies 
the choice f(x) =0. 

The distribution represented by (17) is illus- 
trated in Fig. 1 for several values of ¢. Variations 
in ¢ are considered here to mean variations in R. 
For t=1, the adsorbing material at the end of 
the tube is about 63 percent saturated after the 
given volume has passed. it will be noted that 
only for t>100 do we obtain a close approxima- 











4 rn ie 
' ’ ep 





os 1.0 1.5 , 


Fic. 1. Distribution for k, =0. 
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tion to the discontinuous adsorption predicted 
by the equilibrium theory; that is, only if the 
time of passage of the solution is of the order of 
100 times the “‘half-life’’ of the chemical reaction. 


k, =k, 


For k»=k;, corresponding to an equilibrium 
constant of unity for the reaction, Eq. (11) 
reduces to the linear form 


(18) 


There are no difficulties involved in writing 
down general solutions for this equation, but 


0q/dt+0q/dx+0*q/daxdt=0. 











Fic. 2. Distribution for k, = ky. 


fitting the boundary conditions is a less trivial 
matter. For this reason the Laplace transforma- 
tion was employed.? By the substitution 


q=e-**F(t, x), (19) 
the equation is reduced to 
e°F/daxdt— F=0, (20) 


which leads to a tabulated inverse transform. If 
we let f(s, x) be the transform of F(t, x) with 
respect to the variable /, the transformed problem 
becomes 


sof/dx—f=0; f(s,0)=1/s(s—1), (21) 
which leads to 


f(s, x) =e!*/s(s—1). 


The inverse transform of e*/*/s is Io{2(xt)+} where 


(22) 


2 The notation used here is that of R. V. Churchill, 
Modern Operational Mathematics in Engineering (McGraw- 
Hill Book Company, Inc., New York, 1944). The necessary 
theorems and the appropriate function-transform pairs 
will be found in this text. 
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Io is the Bessel function 


To{2(xt)*} =1-+-axt+ (x?#/[2!)+---. (23) 


The inverse transform of f(s, x) is therefore 


Fit, »)={ e'-TIy{2(x7)*}dr. (24) 


It can be verified without difficulty that the 
solution represented by (19) and (24) satisfies 
all the boundary conditions. This solution can be 
given in series form as 





cs) no) {™ x” 
q=e** 2. > — (25) 
n=0 n+1 m! n! 
~ n t™ )x” 
=1-e'*> | > —}—. (26) 
n=0 ' m=0 m! n! 





Figure 2 includes curves calculated from these 
series for =5 and ¢=10. For larger values of f¢, 
a satisfactory asymptotic solution can be found. 
We write (19) as 


g=i—e=f elt 2(er) "dr. (27) 


Now for large values of the argument, J) becomes 
equal to 





1 
ro exp [2(x7)#]; xr. (28) 


mw (x7 


1 
With this approximation to Jo, g can be written as 


dr. (29) 





1 = [—(xt— ri)?] 
al bane 
2ni Jd, \x ri 
Aside from the factor 7, this integral has the 
form of the probability integral or error function. 
Since this factor is a slowly varying function, we 
will remove it from the integrand. For t<x, the 
integrand has its maximum value at r=x; for 
t>x, the integrand has its maximum value at 
7=t. Using these values for tr when we remove 


7: from the integrand, we obtain 
g=1—3(t/x)*{1—erf (t4—x})}; 
g=1—4{1+erf (x?—#)}; 


x<t, (30) 


x>t, (31) 
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where erf (z) is the error function 


2 z 
erf (8) =— J die (-ila «COD 
0 


T 


which is a tabulated function. For ¢=10, Eqs 
(30) and (31) are nearly correct, as indicated in 
Fig. 2, and will of course be more accurate for 
larger values of t. For very large values of t, Eqs. 
(30) and (31) can be replaced by 


q=3+7% erf (t4}—x!) (33) 


because the variation in (t/x)* is not significant 
in the region of interest. Figure 2 includes the 
curve for t=100; Eqs. (33), (30), and (31) give 
nearly identical results in this range. It will be 
noted that this approximate result (33) satisfies 
the differential equation exactly. It does not 
satisfy the conservation condition exactly, but 
does so to a good approximation. The distribu- 
tion for kx, =k; is similar to that for k,=0, but 
the cut-off is considerably less sharp for the same 
value of ¢. In the equilibrium case, the adsorption 
is discontinuous for the entire range k,<k,;. No 
further new qualitative effects are to be expected 
in this range. 


k, >ky 


For ky>k;, the boundary of the adsorbed 
material is not sharp even under equilibrium con- 
ditions. The differential equation for this range is 


oq q 
Lar 
Ox 
0 1 0g 
+[1+6¢-}———-—=0, (34) 
dx 1+ 6q dt 











Fic. 3. Distribution for K =0.5. 
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where K=k;/ky, B=(1—K)/K. We suggest as 
an approximate solution, valid for large ¢, 


g=}++3 erf (2), 
z= {t/K(1+8q)}*— {(1+6q)x}*. 


(35) 


This solution was suggested by the fact that, if 
the quantity (1+ 8g) were a constant, (35) would 
be a solution of Eq. (34). Using (35) we find that 
the first two terms in (34) are of the order f-! 





of 








Fic. 4. Moving boundary for K = 1.0. 


and x~}, respectively; the last term is of the 
order (xt)~}. The first two terms in (34) combine 
to give a term of the order (x#)—?, which is can- 
celled by a term from the second derivatives. 
Equation (35) therefore satisfies the differential 
equation asymptotically as xt-«. Also, as 
t+, (35) gives the correct equilibrium dis- 
tribution 


(1+69)?=t/Kx; 


Equation (35) satisfies the boundary conditions 
to a good approximation for large values of ¢. It 
is not, however, to be expected that the results 
given by (35) are more than qualitatively correct. 
Figure 3 represents the distribution calculated 
from Eq. (35) for =100, K=0.5. 


Kt<x<t/K. (36) 


CONCLUSION 


The two curves in Fig. 3 do not appear at 
first glance to differ significantly. However, the 
solid curve actually resembles very closely the 
equilibrium curve for K =0.4. This indicates that 
the use of the distribution curve for the deter- 
mination of equilibrium constants may give 
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quite inaccurate results unless the solution is 
passed through the column very slowly. Simi- 
larly, the curve for ¢=100 in Fig. 2 resembles the 
equilibrium curve for a K of about 0.8 rather 
than unity. 


10F 
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Fic. 5. Distribution at x =0 for two components with equal 
equilibrium constants but different rate constants. 


Variations in ¢ can also be considered as vari- 
ations in v at constant R. The above results can 
therefore be used to see how the adsorption 
pattern changes as it progresses down the 
column. For k,=0, we find that the boundary 
remains constant in width as it progresses down 
the column. For k,>0, the boundary becomes 
less sharp as it moves along. Figure 4 illustrates 
this behavior for ky; =k,. The successive curves 
are calculated for t=50, 100, 150; R is assumed 
to be constant. 

It is not probable that much useful information 
can be obtained by studying the rate equations 


WALTER 


for multiple adsorption. For small values of R, 
we would expect simply a blurring of the equi- 
librium patterns. For large values of R, there 
exists of course a possibility of separating two 
solutes which have the same equilibrium con- 
stants but different rate constants. Some in- 
formation on this problem can be obtained by 
investigating the variations in the concentration 
of the adsorbed material at the entrance end of 
the tube. This is a standard problem, so the 
algebraic results will not be repeated here, but a 
curve might be of interest in connection with the 
previous results. Figure 5 represents the variation 
in the concentrations of the adsorbed material 
at /=0. The solutes are assumed to be present 
in the solution in equal concentrations, and the 
two equilibrium constants are assumed to be 
equal. The forward rate constant for solute 1 is 
assumed to be ten times that of solute 2. The 
parameter ¢ is the same as that used previously, 
except that the average value of the forward 
rate constants replaces the forward rate constant 
for the single reaction. For this particular situ- 
ation, we note that if the solution is passed 
through the tube at a rate corresponding to 
t>10, we will have essentially the equilibrium 
distribution, that is, no appreciable separation 
of the two components. However, if the solution 
is passed through the tube at a more rapid rate 
corresponding to = 2, component 1 will be largely 
adsorbed at the beginning of the tube, and an 
appreciable separation will have been obtained. 
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The process of emulsion polymerization is observed to occur in the following manner. An 
aqueous emulsion of monomer is prepared using a soap as the emulsifier. The addition of a 
catalyst starts the polymerization after an initiation period whose length depends on the mean 
radius of the monomer globules (long periods for large mean radius and short periods for small 
mean radius). As the reaction proceeds, the mean radius of the globules diminishes and the 
free soap concentration increases. It is believed that the initiation period is the result of the 
existence of an inhibitor in the monomer. The catalyst reacts with that part of the inhibitor 
which is dissolved in the water phase, so that in order for equilibrium to be restored some of 
the inhibitor must diffuse out of the monomer. Rate curves are derived on this basis using the 
added assumptions: (a) polymerization occurs at the monomer-water interface after the 
inhibitor concentration has reached a low threshold value. (b) radius distribution function of 


monomer has a single maximum. 





A. INTRODUCTION 


ONSIDERABLE theoretical and _ experi- 
mental work has been done on the mecha- 
nism of polymerization in homogeneous systems, 
such as pure liquids or solutions; and it has been 
possible to express the course of such reactions 
with fair approximation by a set of equations 
which represent the different cooperating ele- 
mentary steps of the whole process, such as 
activation, propagation, chain transfer, termina- 
tion, etc. Recently, polymerization in hetero- 
geneous systems, such as emulsions and suspen- 
sions, has become of considerable practical 
importance; and a somewhat detailed study of the 
fundamental background of such processes ap- 
pears, therefore, to be of interest. Only few papers 
(all of them very recent) report experimental 
facts about emulsion and suspension polymeriza- 
tion,! but none of them makes an attempt at a 
qualitative theoretical analysis. 

This article intends to show how, under given 
assumptions, the case of polymerization in emul- 
sion can be treated mathematically. It involves 
such simultaneous processes as diffusion, solu- 
bilization, and chemical reaction, which take 
place in different parts of the system and to- 
gether determine the rate of conversion of the 


1Fryling and Harrington, Ind. Eng. Chem. 36, 114 
(1944); J. R. Vinograd, L. L. Fong, and W. M. Sawyer: 
Abstracts of the papers presented at the ACS Meeting in 
New York City, September 1944, page 3E; W. P. Hohen- 
stein, S. Siggia, and H. Mark, India Rubber World 111, 
173 (1944); 112, 436 (1945). 


monomer into the polymer. It is evident that 
different kinds of assumptions can be made 
regarding the locations of the various elementary 
steps, and it is to be expected that at the present 
incomplete state of our experimental knowledge 
the special assumptions made in this paper will 
require revision and improvements. But it can 
be hoped that the general mathematical pro- 
cedure will, with certain modifications, be useful 
and applicable after more reliable experimental 
data are available. 

Let us, therefore, consider the polymerization 
of styrene in aqueous emulsion. The conditions 
are, in general, the following: About 20 ml 
monostyrene are added to 180 ml water, which 
contains about 1 gram soap (ammonium oleate) 
and with the aid of moderate shaking, a stable 
emulsion of monostyrene in water is produced. 
Previous observations and measurements? have 


TABLE I. Distribution of particle size. (From Clayton’s 
data, reference 3.) 











Soap % % % © % % % % 
concen- 2.0- 2.5- 30- 35- 40- 45- 50- 5.5- 
tration 2.54 3.0n 3.54 404 4.54 S50n S554 6.0 
0.0005M —- —_—_-_ — 14 23 31 18 9 
0.0010 —_-_ — 5 22 34 23 6 6 
0.0015 = 5 22 20 3830 14 9 — 
0.0020 = 18 24 29 18 11 —_ — 
0.0025 6 29 39 13 12 3 —_—_ — 








~ 2 Compare W. Clayton, reference 3; S. Berkman and 
G. Egloff, Emulsions and Foams (Reinhold Publishing 
Company, New York, 1944); and particularly J. W. 
McBain. 
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shown that in such cases a multitude of small 
hydrocarbon droplets are formed, the radii of 
which cover a wide range and exhibit average 
values in the neighborhood of one micron. The 
average particle diameter and the distribution 
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Fic. 1. Variation of mean globule radius of an emulsion as 
a function of emulsifier concentration. 


curve depend upon various parameters, such as 
temperature, mode and vigor of shaking, and 
particularly upon the ratio hydrocarbon (mono- 
styrene) to soap solution (concentration of soap 
and total volume of aqueous phase). Table I 
shows as an example a few rough particle size 
distribution functions*; and Fig. 1 shows how, in 
the case of monostyrene emulsions, the average 
radius of the dispersed particles depends upon 
the soap concentration, while the ratio of hydro- 
carbon phase to aqueous phase is kept constant! 
(1:10). 

These results are commonly explained by 
assuming that those soap molecules which are 
not incorporated in either type of soap micells 
(laminar molecular micells and spherical ionic 
micells), but are monomolecularly distributed in 
the aqueous phase, either as undissociated 
(neutral) acid molecules or as dissociated soap 

* Compare W. Clayton, Theory of Emulsions and Their 


Technical Treatment (The Blakiston Company, Philadel- 
phia, Pennsylvania, 1938), fourth edition, p. 132. 


ELLIOTT W. 


MONTROLL 

and sodium ions, are adsorbed on the surface of 
the hydrocarbon (monostyrene) globules. The 
long (Ciz—C20) hydrocarbon chain of the fatty 
acid (ion or molecule) is soluble in the hydro- 
carbon of the dispersed particles, while the 
COOH, COONa or COO~ groups are highly 
polar and therefore preponderantly attracted by 
the aqueous phase. 

Emulsions of this nature are used to carry out 
polymerization. The reaction is initiated by a 
water soluble substance of peroxide nature called 
the ‘‘catalyst.’’ In most experiments reported in 
the literature, sodium or ammonium persulfate, 
perborate, or percarbonate have been used. Upon 
addition of the ‘‘catalyst”” and continuous agita- 
tion of the system, reaction does not set in at 
once, but one always observes a shorter or longer 
inhibition period during which practically no 
polymer is formed, but some of the ‘‘catalyst’’ 
is used up. This fact is usually explained by 
assuming that certain dissolved impurities in 
the monomer phase or in the water inhibit the 
formation of long chains and must be destroyed 
before polymerization can effectively set in. The 
temperature range in which such reactions are 
carried out is usually between 40° and 80°C. 

Three zones of reaction are possible: 1. in the 
monostyrene phase, 2. in the water phase, and 
3. at the styrene-water interface. 

Since the catalysts used in styrene polymeriza- 
tion are generally inorganic salts insoluble in 
styrene, it is highly improbable that the reaction 
could occur inside the styrene globules. If the 
reaction zone were in the aqueous phase, the rate 
should be independent of total interface surface, 
which does not seem to be the case. When the 
reaction zone is at the water-styrene interface, 
the change in volume of a styrene globule is 
proportional to the surface area of the globule, or 


dV /dt= —4rkr’, 


where & is the proportionality constant. Since 
the sphere volume is given by V=47r?/3, this 
leads to dr/dt=—k, and the radius diminishes 
at a constant rate. This law of disappearance of 
monomer was observed (microscopically) in soap 
solutions of many organic liquids, including 
styrene, by Vinograd, Fong, and Sawyer. In this 
paper we shall consider only those reactions 
which occur at the interface. 
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The rate at which monomer is polymerized in 
a stable emulsion depends on the initial size dis- 
tribution of the suspended globules and on the 
state of agitation of the system. We shall devote 
our attention to initial radius distributions with 
a single maximum, and to reactions which 
proceed while the system is in a state of constant 
agitation. In reactions occurring under these 
conditions one would expect the radius dis- 
tribution curve of monomer spheres to have a 
single maximum at all times. 

In Fig. 2 are plotted experimental curves deter- 
mined by Siggia,! Hohenstein, and Mark, showing 
the rate of polymerization of monostyrene when 
ammonium oleate is used as an emulsifier and 
the reaction is catalyzed by sodium perborate. 
The polymerization was carried out in a shaker 
with 78 shakings per minute. The curves from 
right to left correspond to initial distributions 
having decreasing mean radii. Two interesting 
features of these curves are: 1. The inhibition 
period increases with the initial mean radius, and 
2. The rate curves are essentially parallel in the 
range of intermediate degrees of polymerization. 

Unfortunately the nature and concentration 
of the inhibitor were not determined in these 
experiments. 

Using the properties of monostyrene-water 
emulsions given so far, we shall now discuss a 
polymerization mechanism which qualitatively 
accounts for the curves given in Fig. 2. First‘ a 
process leading to an inhibition period which 
increases with the initial mean radius will be 
suggested, and then by employing initial dis- 
tributions of the type given in Table I we shall 
describe how this process would yield rate curves 
similar to those found in Fig. 2. 

If it is assumed that the initiation period is due 
to the presence of an inhibitor in the absence of 
a catalyst, this inhibitor is distributed between 
the water and styrene so that the system is in 
equilibrium. As soon as the catalyst (which in 
general is much more soluble in water than in 
styrene) is added, it reacts with the inhibitor in 
the water, destroys its equilibrium between the 
two phases, and causes inhibitor to diffuse from 
the styrene into the water. If there exists a 
critical concentration of inhibitor in the styrene 


* This process is similar to that given by Mark, Hohen- 
stein, and Siggia. 
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below which it polymerizes, it is clear that poly- 
merization would start in a styrene globule as 
soon as enough inhibitor diffuses from it to 
reduce its inhibitor concentration to a value less 
than the critical. Since diffusion from a large 
globule is slower (as will be shown in the mathe- 
matical section of this paper) than from a small 
globule, the small globules in an emulsion would 
start to polymerize before the larger ones, and 
emulsions of small mean globule radius would 
polymerize faster than those of large mean 
globule radius provided that both have the same 
initial inhibitor concentration. This is in agree- 
ment with the rate curves in Fig. 2. 

If initially a styrene-water emulsion has a 
sphere size distribution similar to one of those in 
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Fic. 2. Rate of polymerization of monostyrene for 
various initial mean globule radii when ammonium oleate 
is used as an emulsifier and sodium perborate as a catalyst. 














Table I, according to our model, upon addition 
of catalyst the smaller globules start to poly- 
merize first. Since there are initially relatively 
few small globules, the extent of polymerization 
changes very slowly at first. As time goes on the 
larger globules become available for polymeriza- 
tion as the result of two effects. First, more 
inhibitor diffuses out of them leaving them with 
inhibitor concentrations lower than the critical. 
Secondly, as the small globules polymerize, they 
give up their soap, thus increasing the emulsifier 
concentration and reducing the sizes of the 
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larger globules, causing them to polymerize 
sooner than they would otherwise. 

One would expect both of these effects to con- 
tribute to an increasing polymerization rate 
until practically all the globules are able to start 
polymerizing, after which time there could be no 
increases in the rate. This process is in agreement 
with the individual rate curves in Fig. 2. 

In the next part of this paper we shall make a 
detailed mathematical study of the processes 
described here. It will be divided into six 
sections: 


I. Inhibition time. 
II. Variation in mean radii with variation in soap con- 
centration. 
III. General reaction equation. 
IV. Globule radius distribution functions. 
V. Solution of reaction equation. 
VI. Discussion of numerical results. 


The effect of swelling of polymer by residual 
monomer on reaction rates is neglected in the 
following discussions. 


B. MATHEMATICAL SECTION 
I. Inhibition Time 


One of the most striking properties of the 
polymerization rate curves in Fig. 2 is the 
initiation period which depends on the initial 
mean radius of the spheres in the dispersed phase. 
We have already suggested that this period in a 
particular globule is related to the time required 
for a certain amount of inhibitor to diffuse out 
of the globule of monomer. In order to apply this 
conjecture to a mathematical theory of the rate 
curves, we must give a definite description of the 
diffusion process. In particular we must derive 
a relation between the radius of a suspended 
monomer sphere and the time required for the 
sphere to achieve a state from which polymeriza- 
tion can start. 

Let us suppose that at time t=0 a sphere of 
radius R of one liquid, Lo, is suspended in another 
liquid, L., and that at equilibrium Co is the con- 
centration of a solute S in Ly and c¢, is the con- 
centration of S in L,. If to the L.—Lo system a 
reagent is added which reacts with the solute 
dissolved in L. but not that dissolved in Lo, 
equilibrium is destroyed, and solute must diffuse 
into L, until equilibrium is reestablished. Under 


the assumption that catalyst reacts with the 
solute S as soon as it diffuses into the water phase 
(i.e., ce=0 if 20) we shall calculate the concen- 
tration of S at various points in the sphere Lo 
as a function of time. 

According to Fick’s law, the concentration 
c(r, t) of S at a point a distance 7 from the center 
of the sphere at time ¢ is given by 


D(6?/dr?) (rc) = (d/dt) (rc), 


where D is the diffusion constant of S in Lo. If 
at time t=0 the sphere has a uniform solute con- 
centration ¢ , we have the initial conditions when 


t=0, r<R, (rc)=reo, 


and if there is infinite permeability at the L»)—L,. 


interface, we have the boundary conditions: 
when ¢>0, r=R, (rc)=0; r=0, (rc)=0. 


One can easily show that the solution of the dif- 
fusion equation under these conditions? is 


—]1 n+1 


— sin (nrr/R) 





c(r, t) =(2Reo/zr) > 


n=1 n 
Xexp — {n?x*Dt/R?}. 


If C(t) is the average concentration at time ¢, 


C(t) = (3/4aR) f Anr’c(r, t)dr 


= (6€0/7) > (1/n?) exp — {n?x*Dt/R?}, 


n=1 


and as ¢ increases only the first term becomes 
significant : 


C(t)/com(6/m) exp —2’Dt/R?. 


So, the time required for the average concen- 
tration of a sphere of radius R to become C(é) is 


t~(R?*/D*) log (6¢0/rC(t)). (1) 


In our theory of emulsion polymerization we 
shall assume that there is a critical average con- 
centration c* which a globule must have before 
it polymerizes. Thus, a globule of radius R will 
not start to polymerize until the time ¢ given by 
(1) when C(t)=c*. Stated in another way, at 


5W. Byerly, Fourier Series (Ginn and Company, New 
York, 1893), p. 116. . 


nr -« a 
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time ¢ all spheres of radius less than 
R~#/{(1/mD) log (6co/mc*)}*, (1a) 


will be polymerizing. This expression can be 
written in a simpler form by introducing the 
time, 7, required for a sphere of unit radius to 
start polymerizing. In terms of 7, we have 


R™(t/r)?. (2) 


There are many objections that could be 
raised in the manner in which (2) was obtained, 
but even under more exact treatment one finds 
that R is still roughly proportional to #}. 


II. Change in Sphere Diameter with Emul- 
sifier Concentration 


It has been observed experimentally that the 
mean radius of the suspended spheres in an 
emulsion decreases with increasing emulsifier 
concentration. In Fig. 1 we have plotted the 
mean sphere radius of a styrene emulsion in 
water as a function of emulsifier (ammonium 
oleate) concentration. The curve indicates that 
(except possibly in the range of very small 
radii) the mean radius varies linearly with the 
concentration, and can therefore be represented 
by the equation 


(r —1ro)/(C—Co) = (ro— 171) / (Co— C1), (3) 


where (79, Co) and (ri, ¢1) are two points on the 
linear portion of the curve. In our analysis it 
will be assumed that the mean radius, 7, is a 
linear function of the concentration over the 
entire concentration range. 

From (3) it is possible to obtain a relation 
between mean radius and extent of reaction. We 
define the extent of reaction as 


a=[V(0)—V()]/V(), (4) 


where V(t) is the total volume of suspended 
spheres of monomer at time ¢ (we assume that 
the spheres of monomer have a constant density 
independent of radius). Another convenient 
variable is 


B=1-a. (4a) 


Suppose Co is the emulsifier concentration cor- 
responding to the initial mean radius 79 at the 
beginning of the reaction. Since in the reactions 
in which we are interested emulsifier is adsorbed 
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on the monomer globule surfaces and is given up 
by the monomer as it is polymerized, the free 
soap concentration increases linearly with the 
loss in total surface of the monomer globules, 
thus 

c=cot+c’(s(0) —s(t)), 


where c’ is a constant and s(t) is the total surface 
of all the globules at time ¢. If r1;=0, ci is the 
emulsifier concentration when the reaction is 
completed (i.e., when s(¢) =0) or 


C1=Cot+c's(0). 


Introducing these concentration relations into 
(3) we have 


(r —10)/(s(0) —s(t))c’ = —10/c’s(0), 
r=ros(t)/s(0). (4b) 


Unfortunately the ratio s(t)/s(0) is not a con- 
venient quantity for later calculations; however, 
if the ratio V(t)/V(0)=6 could be used to 
approximate it, the equation 


r—roB (5) 


would be useful in developing our reaction-rate 
theory. 
When the reaction begins, :=0 and 


s(t)/s(0) = V(t)/V(0) =1, 
and when the reaction is completed 
s(t) = V(t) =0. 


So (5) is exact in the two extremes 8=1 and 
8=0. For the binomial and normal globule 
radius distributions that we shall use later, the 
agreement between the volume ratio and surface 
ratio turns out to be as good as necessary for the 
calculations in this paper. The largest percentage 
disagreement occurs in the range of small mean 
radii, where the assumption of a linear relation 
between soap concentration and mean globule 
radius is no longer valid anyhow. On the basis of 
these approximations our rate curves should be 
least accurate at the end of the polymerization 
process. 

The error involved in the application of (5) 
is probably largest in the case of a sharp radius 
distribution. In fact, with a uniform radius dis- 
tribution (all spheres with the same radius) we 
have 
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Fic. 3. Sphere size distributions for various emulsifier concentrations as taken 
from Clayton’s data. 


s(t)/s(0) =N()r°/N(O0)r0’, 
and 
B= V(t)/V(0) =r®N(t)/ro®N(0), 


where NV(t) is the total number of globules (all 
with radius 7) at time ¢. Thus, in this case, com- 
bining the above two equations with (4b) 


r=7r 6). (5a) 


In broader distributions the exponent on 8 is 
probable between 3 and 1. 

One immediate result of (5) and (5a) is that 
they allow one to express the mean radius of the 
monomer-polymer mixture as a function of the 
extent of the reaction, a. If we let 


m= Mean monomer radius, 
’,»=mean polymer radius, 


the mean radius of hydrocarbon spheres that 
would be observed experimentally is (using (5)) 


(rw =%m(1—a)+ar, 


=r79(1—a)?+ar,. 


(6) 


A similar result can be obtained in the case of a 
homogeneous distribution from (5a), 


(r)w=ro(1—a)!-tary. (6a) 


If for all extents of reaction the mean radius 
of polymer globules is the same as their final 
mean radius after polymerization is completed, 
(6) and (6a) allow one to calculate (r)w from a 


knowledge of the initial and final mean radius 
of the globules. Early in the reaction this as- 
sumption is probably in error due to the swelling 
of polymer globules by monomer; however, 7, has 
a as a coefficient so that the mean radius is less 
sensitive to errors in 7, in this range. 


III. General Reaction Equation 


We shall now derive a differential equation for 
the extent of reaction under the assumptions that 


a. Reaction takes place at boundary of suspended 
spheres, 

b. Reaction rate is proportional to total surface area of 
suspended spheres, 

c. No sphere begins to react until its average inhibitor 
concentration reaches a critical value, say c*, and 
after that the reaction rate of a particular sphere is 
independent of the inhibitor concentration of that 
sphere. 


As was shown in section I, at a time # all 
spheres of radius less than R are polymerizing if 
R is given by R=(t/r)!, where 7 is the time 
required for a sphere of unit radius to start 
polymerizing and is related to c* by (1a) and (2). 


TABLE II. Effect of emulsifier concentration on mean 
radius and standard deviation of sphere radius 
distribution function. 








Soap conc., 10 M 5 10 15 20 25 
Mean radius, 7, 
(microns) 4.67 4.36 4.01 3.65 3.27 


Standard deviation, ¢ 0.587 0.600 0.672 0.620 0.587 
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If G(r, t)dr is defined to be the number of 
monomer spheres of radius between r and r+dr 
at time ¢, assumptions a—c and Eq. (2) imply that 
the rate at which the total volume of monomer 
is suspension changes at time ¢ is 


(t/r)* 
dV(t)/dt=—k f 4nrG(r, t)dr, (7) 
0 


where k is the rate constant for the reaction. 
Now suppose JV(#) is the total number of spheres 
in suspension at time ¢. Then the fraction of 
spheres with radii between r and r+dr at time 
tis 

g(r, t)dr=[G(r, t)/N(t) Jdr. 
Clearly, 


N()= f “Gt, t)dr, 


V(t) = f ” (4/3) arG(r, t)dr 
=Ny f (4/3)"40, t)dr. 
0 
Since from (4) V(t) = V(0)8, 
~(d log B/dt) J “g(r, dr 
= 3k fe, t)dr, (8) 


and, given a sphere size distribution g(r, t), one 
can calculate the extent of reaction at a given 
time from this equation. Of course g(r, ¢) is a 
function of 8, so (8) is really a differential- 
integral equation. For calculating purposes it is 
convenient to introduce two new parameters, T 
and \ defined by 


T=3kt, \=3kr, (9) 
which change (8) into 


~(1/8)(d8/aT) f g(r, t)dr 
; (t/r)4 
~ f Pg(r, fdr. (10) 


Equation (10) can be considered to be the funda- 
mental reaction rate equation for emulsion 
polymerization. 
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IV. Sphere Size Distribution 


As was mentioned in the general discussion, 
the sphere radius distribution of a stable emul- 
sion can be expected to have a single maximum 
provided the emulsion is well stirred. This was 
exemplified in Clayton’s data in Table I. The 
distribution functions of Table I are plotted in 
Fig. 3, as are their approximations by normal 
curves; these normal curves seem to be good 
approximations to Clayton’s data. 

The decrease in mean radius with increasing 
emulsifier concentration is clear from Figs. 1 
and 3. In Table II we give the mean radius and 
standard deviation of globule radii from the 
mean as obtained from Clayton’s data. This 
table indicates that to within the accuracy ex- 
pected in this work the mean radius, #, is a linear 
function of the soap concentration and the 
standard deviation is almost independent of the 
soap concentration. In the distributions given 
in Table I, the mean value of o is 0.613y. 

On the basis of the data given in Tables I and 
II we shall for our numerical work assume that 
initially the radii of monomer spheres are dis- 
tributed normally or according to binomial ap- 
proximations of normal distributions. Further- 
more, we shall assume that the standard devia- 
tions from 7 of these distributions are constants 
in a given emulsion, independent of the soap 
concentration. 

If g(r, t)dr is the fraction of monomer spheres 
with radii between r and r+dr at time ¢, a 
normal distribution of sphere radii can be ex- 
pressed analytically by 


g(r, t) =[1/o(2m)*] exp —(r—7)?/20?, (11) 


# being a function of the time and according to 
(5) expressed by 


Fo~ro(1—a) = 708, 


where a is the extent of reaction a time ¢ after 
catalyst has been added to a monostyrene emul- 
sion with an initial mean sphere radius ro. The 
simplest binomial approximation to a normal 
distribution is a rectangular distribution with 
the same mean radius and standard deviation. 
Analytically this can be written: 


Ken 0 if |r—#| >o3! 
B's 1/(2034) if |r—#| <o34. (12) 
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A more complicated but more accurate binomial 
approximation to a normal function is the three 
rectangle approximation which can be repre- 
sented analytically by 


() if (r—#) > 36/2 

1(1/46) if 36/2>(r—7) >8/2 

g(r, )=42(1/48) if 6/2>(r—*A) >—6/2 ~— (13) 
1(1/48) if —8/2>(r—7) >36/2 

0 if —36/2>(r—7), 





where 


§=40(3/7)}. 


By taking more and more rectangles one can 
obtain increasingly more accurate binomial 
approximations. 


V. Solution of the Reaction Equation 


In this section we shall solve the reaction 
equation for several initial globule distributions. 
The physical nature of the solutions will be dis- 
cussed in Section VI. Let us first consider the 
case of no initiation time. This might be ap- 
plicable to situations where inhibitor is initially 
on the surface of the hydrocarbon globules but 
not in their interiors. 

Whether there is an inhibitor or not in the 
globules, one must, in order to solve the reaction 
Eq. (10) evaluate 


(73) = f “Pele, t)dr. 


This we shall do in three cases. 
(a) Rectangular Distribution. Application of 
(12) yields . 


? +034 
(A) = A(1/203')dr 
?—o34 


= (1/803*) { (7+-03!)4— (7 —03!)} 
= 7(7+30°). 


(b) Binomial Distribution. (6 =20(3/7)*). Ap- 
plication of (13) yields 


y= i) rye? 8 /48)dr-+2 f ree 8/48)dr 


7 +6 /2 7-5/2 


#-6/2 
+ f (r?/46)dr 
7-36 /2 


= 7(7+30"). 
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(c) Normal Distribution. 


()= f C/o(28)" exp —(r—7)?/20°%dr 
“ = 7(7+36°). 


It is to be noted that in all three cases (and 
indeed for any binomial distribution) 


(3) = 7(7 +30). (14) 


With no initiation period we must use 
(")= f rg(r, dr=F+o% (15) 
a ; 


With (14) and (15) one can solve (10). 


1. No Initiation Period 


When there is no initiation period, (10) be- 
comes, upon application of (14) and (15): 


— (7/8) (d8/dT) (7 +30") = (0? +7"), 


a result valid for all binomial distributions. 
Since *= ro, 


eee 
ital nih argumiee 
B+ (a/ro)* 


Since at time ¢=0 (i.e., 7=0) no polymerization 
has occurred, 8=1, and the solution of (16) is 


| =ar. (16) 


T =aro+2o(tan ro/o—tan— Bro/a). (17) 


This equation is plotted in Figs. 4 and 5 for 
several values of 79 and o. 


2. Homogeneous Distribution with Inhibitor 


In the presence of an inhibitor the initial 
radius distribution easiest to treat is a homo- 
geneous one. From the discussion in Section III, 
8=1 in this case until the time corresponding to 


T= To=Aro’,: 


where 7» is the initial radius of all the hydro- 
carbon globules. After the time 7», 6 is the 
solution of the differential equation obtained by 
introducing (5a) into (10) and choosing 


1 if r=r6? 
0 otherwise: 


— (1/8) (dB/dT) (roB*)* = (70884)? 


g(r, t) = 
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Fic. 4. Theoretical polymerization rate curves in the 
absence of inhibitor. Standard deviation of sphere radii 
from mean 0.6p. 


or 


—r dB/B'=dT. 


The solution of the equation is (employing the 
boundary condition B=1 when T=T>=Ar¢*) 


1-—a=8=[1—(1/2ro)(T—Are) ?. (17a) 


This equation is plotted in Fig. 6 for several 
values of 79 when A = 1. It is to be noticed that the 
reaction starts suddenly at the end of the 
initiation time instead of gradually as one ob- 
serves experimentally. By assuming a sphere 
size distribution of the type discussed in Section 
IV, one obtains theoretically the observed 
gradual increase in reaction rate at early times. 


3. Rectangular Distribution with Inhibitor 


The number of globules being polymerized at 
a given time in this case depends on the time. 
Since there are initially no globules of radius 
less than (79—03'), there can be no polymeriza- 
tion at times less than those corresponding to 
To/X=(ro—03')?. For values of T/X between 
T)/X and T1/A= (1708103)? only part of the 
globules polymerize. Finally for times >7,/) all 
globules have a sufficiently low inhibitor con- 
centration for polymerization to take place. 
During each of these three time intervals the dif- 
ferential Eq. (10), 


(T/r)t 
ro(dB/dT) (746?-+ 302) = — f Pg(r)dr, 


0 
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Fic. 5. Theoretical polymerization rate curves in the 
absence of inhibitor. Standard deviation of sphere radii 
from mean 0.3. 


has a different form. Applying (12), 


(a) When TST», d8/dT=0; (18a) 
(b) When TMo<TS Ti, 





1 (T /r)t 
ro(dB/dT) (r°8? +30?) = — J rdr, 
203% Bro—o(3)4 
or 


ro(d8/dT) (10°68? +307) 





1 
= ———[(T/d)!—(ro8—03*)*]; (18b) 
603 


a3} 
(c) When T>T7;, 
ro(dB/dT) (10°68? +30") = —(r0°8’+0*). (18c) 


Thus in each time interval the form of the solu- 
tion of (18) is different. These solutions are to 
be combined through the boundary conditions 
at the beginning and end of each interval. 
Interval (a). T/X<To/d. The solution of (18a) 
is 8=constant. To fit the boundary condition of 
no polymerization at t=0, the constant must be 


1, so for 
T/XX£To/r, B=1. (19a) 


Interval (b). To < T<T;. In this interval the 
solution of the differential Eq. (18b) is rather 
complicated. We shall assume that 


B(T)=1+> a,(T—T>)"/n!. (19b) 
n=1 
According to Taylor’s theorem, 


dn=d"B(T»)/dT™. 








If we let ws=7o0(7e?+30?) and remember that 
8=1 when T=Tpo, we find by successive differen- 
tiation of (18b) that 


a;=d8/dT]7_7,=0, 
a= To'/[4odus(3d)*], 
a3>— { (To?/4A3) —de2ro(To/d) } /2ou33?. 
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Fic. 6. Theoretical variation of degree of 
polymerization, a, with time. 


For reasonable values of 79, J, A, and o the main 
contribution (over 90 percent) of a3 comes from 
the term proportional to a2, and in general, the 
dominant term in @, is proportional to a,_1, so 
that 

On—~An—17 0(T0/d)/ (20334). 


Application of this approximation to (19a) yields 
B—1~— {opus(3dr)*/re?T 93} 
X {exp (T— To) (r0T 0/20Au334) —1 

i (T- To)roT o/2opn3d3'+0(T— T»)?} ° 
The correction terms of O(T—T>)* can be 
written as a power series in (T—TJp) (starting 
with (T—T>)*). Numerically they are small 
compared with the first three terms in the above 


equation; however, for completeness we shall 
include several of them in the final expression :. 


—a=B—1= — {ops(3d)?/r°T 03} 
X {exp (T— To) (roT0/20dAu33*) — 1 
~(7~7Z) T oro/2odp33* 
+ Tore?(T — T0)?/14407 u 3"? 
+[27?To?(3Todr)4/od2 3? —1 |ro? 


X (T—To)*/11520°r?us?+---}. (20) 
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This expression is valid until JT becomes suf- 
ficiently large so that 

T =)[ro8 +o0(3)* }. (21) 
For larger values of T, (18b) is no longer valid, 
and (20) does not apply. If we define 7; and f, 
as the values of T and 8 for which (20) and (21) 
are satisfied simultaneously, we have as boundary 
conditions for (18c) T=7, when B=A;. 

Interval (c). T>Ty. In this interval one can 
easily obtain a closed relation between T and 6 
by solving (18c) with the above boundary con- 
ditions. Thus 


(Ti1—T) =ro(8—B:) 
+2¢(tan Bro/o—tan-! Biro/c). (22) 


The complete solution of (10) for all times is the 
combination of the solutions in intervals a, }, 
and c. This solution is plotted in Figs. 7-8 for 
several values of 79, o, and X. It will be discussed 
in detail in Section VI. 


4. Binomial Distribution with Inhibitor 


In the binomial distribution (11) there is 
initially no globule of radius less than ro9— (35/2) 
(where &6=40(3/7)!) so there cannot be any 
polymerization until times corresponding to 
To=X(ro— 35/2)”. With this distribution we 
must consider separately five time intervals in 
each of which (10) takes on a different form. 
Whenever g(r, ¢) changes its form a new interval 
begins. If these intervals begin at times cor- 
responding to To, 71, T2, ... we shall let the 
values of 6 at these times be 1, 81, Bo, .... The 
five required differential equations with their 
boundary conditions are: 

(a) When T<To=X(ro—36/2)?, 


dp/dt=0 (23a) 


with the boundary condition B=1 when T=0. 
(b) When Bro—36/2<(T/d)!<Bro— 5/2, 


ro(d8/dT) (70°78 +-30”) 
= (—1/126) {(T/d)?— (708 —36/2)*} (23b) 


with boundary condition B=1 when T=T» 
= Aro 36/2)?. 

(c) When Bro— 6/2 < (T/A)? <ro8+56/2, 
ro(70°?8?+307)dB/dT 


= — (1/126) {2(T/A)!— (Bro—8/2)* 
— (Bro—36/2)*} (23c) 











wi 


ro 


suf- 


(21) 
alid, 
id Bi 

(21) 
dary 


can 
nd B 


con- 


(22) 
the 
1, b, 
for 
ssed 


> is 
5/2) 
any 
to 
we 


rm. 
‘val 
-or- 
the 
The 


1eir 


3a) 


3b) 
Te 


3c) 








with the boundary conditions 8=8,; when T=T7, 
= (7081 — 5/2). 
(d) When 6ro+36< (T/A)! <Bro+36/2, 


ro(B’ro?+30*)dB/dT 
= — (1/128) {(T/d)?+ (ro8 +38)? 
— (ro8 — 36)* — (ro8—35/2)*} (23d) 


with the boundary condition B=82 when T=T>» 
= (Boro +35). 
(e) When (T/A)!>6r,.+36/2, 


ro(r0°B’+30°) (dB/dT) = —(ro°h’+o0°) (23e) 


with the boundary condition 8=83; when T=T3 
= )(B3r0+36/2)?. 

In the interval (a), as in the case of a rectan- 
gular distribution, 8=1 when T<7>. In interval 
(b) we proceed in the same manner as in the 
rectangular distribution and obtain: 


To(T— To) 
exp 
46d(1r70° +30?) 





25us(A/To)? 
re°To 
T.(T — To) 
45X(r02-+302) 


—a=p-1=— 








+0O(T—T>)* . (24a) 


In interval (c): 
B=6it+a,(T— T;)+a2(T— T,)?/2!+- es (24a) 


where 





1 — — 36(T;/r)' +8? 


a= —— 
re BYr+3o" 


, etc., (24b 
1275 ) ( 


each a; being obtained by taking successive 
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Fic. 7. Theoretical polymerization rate curves with 
inhibitor. 
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derivatives of (23c). In interval (d): 
B=B2+a,(T— T2)+a2(T— T2)?/2!+ did (24c) 


where 





1 /3(T2/d) —58(T2/d)'+38 
( ). etc. 


a= —-— 
4ro ro°Bo? +30? 





| =f 
°o 


e 
merpUurhr uansN & & 
PTET rr dt 


——- 
——S 
—— 
——_— 
—— 
—o 
aa 
—-~ 














Fic. 8. Theoretical polymerization rate curves at con- 
stant initial inhibitor concentration, standard deviation, 
and mean radius. 


In interval (e): 
(T—Ts3) =ro(B—Bs)+20(tan Bro/o 
| —tan-! Byro/a). (24d) 


The complete curve of a as a function of time 
as obtained by combining the above equations 
is plotted in Fig. 7. 

In the case of a normal distribution Eq. (10) 
cannot be solved directly by quadratures, but 
any of the well known methods of numerical 
integration of first-order differential equations 
can be applied. The results are similar to those 
we have obtained from binomial distributions. 


VI. Discussion of Numerical Results 


Equation (17) of the last section gives the 
extent of reaction, a, as a function of time when 
(a) no inhibitor is present, and, (b) the initial 
monomer distribution is binomial. From (17) it 
is clear that the dependence of a@ on time is 
determined by two parameters, 7o, the initial 
mean radius of the monomer globules, and ¢, the 
standard deviation of their radius distribution. 
Figures 4 and 5 show that under conditions (a) 
and (b) polymerization starts immediately, and 
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Fic. 9. Theoretical variation with time of mean 
radius of monomer-polymer mixture. 


of all systems with the same initial mean radii, 
those of the smaller standard deviation (that is, 
those with the sharper radius distributions) 
polymerize faster. Of all those of the same 
standard deviation those with the smaller mean 
radii polymerize faster. 

Equations have been derived for the progress 
of the polymerization reaction when inhibitor is 
present for two initial sphere radius distribu- 
tions: (a) rectangular distribution, and (6) bi- 
nomial distribution. In case (a), @ is given as a 
function of time in Eqs. (19), (21), and (22) while 
the results for case (b) are given in (24a)—(24d). 
Besides the parameters 79 and o, one more 
parameter A, a measure of the time required for 
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inhibitor to diffuse out of a monomer sphere of 
unit radius, is needed to completely describe the 
variation of a with time. @ is plotted in Fig. 7 
as a function of time for several initial mean 
radii when A\=1 and o=0.6. It is to be noted 
that, except for very early times, the rate curves 
of distributions (a) and (b) corresponding to the 
same initial mean radii are very similar. Indeed, 
the detailed shape of the initial distribution is 
not nearly as important as its standard deviation 
and mean radius in determining the character of 


‘the extent of reaction curves. The curves in Fig. 


7 resemble closely the experimental curves of 
Fig. 2 in the manner in which they exhibit the 
increase in the initiation period that occurs when 
the initial mean radius is increased. In Fig. 8 the 
extent of reaction is, plotted for several values 
of \ while o and 7p are kept fixed. The larger the 
value of \ the longer the inhibition period. 

In Fig. 9 the mean radius of the polymer- 
monomer mixture is plotted as a function of time 
using Eq. (6) and the curves corresponding to the 
binomial distribution in Fig. 7. 

In conclusion the author wishes to thank 
Professor H. Mark and Dr. W. P. Hohenstein, 
who have been doing experimental work on 
emulsion polymerization, for suggesting this 
problem and for their helpful advice during its 
investigation. 
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HIS section will accept reports of new work, provided 
these are terse and contain few figures, and especially 
few halftone cuts. The Editorial Board will not hold itself 
responsible for opinions expressed by the correspondents. 
Contributions to this section should not exceed 600 words in 
length and must reach the office of the Managing Editor not 
later than the 15th of the month preceding that of the issue in 
which the letter is to appear. No proof will be sent to the 
authors. The usual publication charge ($3.00 per page) will 
not be made and no reprints will be furnished free. 





Molecular Constants and Chemical Theories 


W. F. Luper 
Northeastern University, Boston, Massachusetts 
July 9, 1945 


T seems obvious that the interesting controversy be- 
tween R. Samuel and C. W. Wheland! cannot be settled 
by appeals to wave mechanics. But citation of a few experi- 
mental facts will show that the extreme position for and 
against the octet theory cannot be maintained by either. 
Samuel’s reaction against the octet theory is justified, 
but cannot be supported by discarding the coordinate bond. 
Wheland’s reaction against abandoning the coordinate 
bond is justified, but cannot be supported by maintaining 
the octet theory. The “‘octet theory” must be abandoned, 
but the coordinate bond cannot be given up. 

The existence of the chemistry of additiori compounds 
(not mentioned by either Samuel or Wheland) requires 
the concept of the coordinate bond, e.g., in the reaction 
between BCI; and NH3. On the other hand, the existence 
of a large number of “‘violations’’ of the octet rule requires 
its abandonment. Covalent compounds such as PFs, SFs, 
IF;, and OsFs show that the octet rule is rigorously valid 
only for the second period of the periodic table from Li to 
Ne. This is admitted by Dr. Wheland in a footnote. 


1J. Chem. Phys. 13, 239 and 251 (1945). 





Solute Type and Viscosity of Dilute 
Solutions of High Polymers 


RUDOLPH SPEISER AND R. T. WHITTENBERGER 


Eastern Regional Research Laboratory, 
Philadelphia 18, Pennsylvania} 


July 2, 1945 


RITH®? has recently considered the theoretical sig- 

nificance of the slope of 9.»/c—c curves for solutions 
of high polymers. By statistical and quasithermodynamical 
analysis, she has shown that the slope is related to the 
energy of polymer-solvent interaction, the slope being 
greater for a good solvent or solvent mixture than for a 
poor one.2* Frith has also shown that at low concentra- 
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tions her equation for the variation of msp/c with ¢ is 
similar to Huggins’.® 

In Huggins’ equation the slope, k’[n }, is proportional 
to the square of the intrinsic viscosity, and the constant, 
k’, is characteristic, and independent*® of the intrinsic 
viscosity for a given polymer-solvent system. This in- 
variance of k’ makes it a more useful quantity for discus- 
sion than the slope. For this reason we have combined 
Frith’s and Huggins’ equations to obtain the fellowing 
relation among the constants: 


k’ = B'/(n P—Bw/[n PRT =A —B"'w/T. (1) 


where A and B” are independent of [ny]. 

The second term accounts for the variation of 9.»/c with 
concentration due to the energy of interaction, w, between 
the polymer chain segment and solvent molecules, whereas 
the first term accounts for variations with concentration 
due to other causes, such as differing internal strueture. 
w is negative? when a tendency exists for solvation of the 
polymer chain segment. Thus, if a solvent is made in- 
ferior by addition of a non-solvent, which leads to a less 
negative w, k’ decreases; if a better solvent is used, w has 
a more negative value and k’ increases. Likewise, if the 
energy of interaction is constant and only the internal 
structure is changed, k’ varies with A. 

Variations of k’ can be divided into two cases: (1) the 
interaction energy, w, remains constant and A varies; 
(2) A remains constant and w varies by changing (a) the 
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FiG. 1. nsp/c —c curves for potato amylose and amylopectin fractions 
in 2 percent KOH. Solid lines amylose. Dashed lines amylopectin. 
Amylose and amylopectin were separated from whole starch with 
nitrobenzene. (See reference 8.) Amylose components (and amylopectin 
components) were fractionally precipitated by isopropyl! alcohol from 
0.5 percent KOH. 
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TABLE I.—X’ values calculated from sp/c —c curves for 
pectin nitrate-acetone solutions. 











Treatment of Degree? 
pectin sample of methyl 
before nitration esterification k’ N 
Untreated 0.75¢ 0.87 +.02¢ 9.4% 
Acid deesterified* 0.24¢ 0.19 +.02¢ 9.4 








@For method, see C. H. Hills, Ji W. White, Jr., and G. L. Baker, 
Proc. Inst. Food Tech. page 47 (1942). 

’ = [COOCHs] /( [COOH] + [COOCH3]); C. H. Hills and R. 
Speiser, Science, in press. 

¢ Degree of methyl esterification before nitration. Nitration results 
in approximately 30 — methyl deesterification. 

4 Average for five fractions of different intrinsic viscosity. 


solvent or (b) the solute type. Case 2(b) is often trivial 
but becomes important when one studies. the change in 
polymer-solvent interaction as the nature of the surface 
of the molecule is changed. 

Numerous examples of Case 2(a) may be found in the 
literature.245 A possible example of Case 1 was the ob- 
servation of Bartovics and Mark® on polystyrene that k’ 
changes with the temperature of polymerization. This 
variation in k’ was attributed to a change in the internal 
structure of the polystyrene as a result of a change in the 
mechanism of polymerization with the temperature. 

We have obtained viscosity data which provide good 
examples of Cases 1 and 2(b). Figure 1 represents ysp/c—c 
curves for fractionated potato amylose and amylopectin 
(Case 1). It is evident that k’ (computed from Huggins’ 
equation) for amylose (linear chain) and amylopectin 
(branched chain) fractions, respectively, is independent of 
the intrinsic viscosity. However, k’=0.58+0.02 for amy- 
lose? is considerably less than k’=1.49+0.10 for amylo- 
pectin. These two substances have a_ well-recognized 
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difference in molecular structure. Since both amylose and 
amylopectin are composed of glucose residues, w for amy- 
lose and amylopectin must be nearly identical, and the 
difference in k’ must be due to a difference in A for the two 
molecules. These results suggest that this method may 
provide a means of detecting branching in high polymers. 

Case 2(b) is illustrated by pectin nitrate solutions in 
acetone for which Table I lists average values of k’ for 
five fractions. Two samples were used, their degree of 
methyl esterification before nitration being 0.75 and 0.24. 
The ’s for the pectin nitrate fractions obtained from ez 
sample are independent of the intrinsic viscosity, but 
for “high ester’ pectin nitrate is greater than that for 
“low ester” pectin nitrate. In this example, the solt 
molecule is principally altered through replacement of 
methyl ester groups by carboxyls, thus changing polymer- 
solvent interaction, with consequent variation of k’. 
Because methyl ester groups are more compatible with 
acetone than carboxyl groups, k’ for pectin nitrate made 
from the untreated pectin is greater than k’ for the acid 
deesterified pectin nitrate. 


1 One of the laboratories of the Bureau of Agricultural and Industrial 
Chemistry, Agricultural Research Administration, United States De- 
partment of Agriculture. 

2E. M. Frith, Trans. Faraday Soc. 41, 17, 90 (1945). 

3E. M. Frith and R. F. Tuckett, Nature 155, 164 (1945). 

( — A. Bartovics, and H. Mark, J. Am. Chem. Soc. 64, 1557 
194 

5M. L. Huggins, J. Am. Chem. Soc. 64, 2716 (1942). 

6 A. Bartovics and H. Mark, J. Am. Chem. Soc. 65, 1901 (1943). 

7 Data of J. F. Foster and R. M. Hixon, J. Am. Chem. Soc. 65, 618 
(1943), yield k’ =0.62, 0.67, 0.57, and 0.65 for whole potato amylose, 
whole tapioca amylose, whole lily bulb amylose, and whole corn amylose, 
respectively. This agreement with our data for amylose is noteworthy, 
since Foster and Hixon’s solvent was anhydrous ethylene diamine, 
whereas we used 2 percent KOH. 

8 R. L. Whistler and G. E. Hilbert, ‘‘Theory of the Fractionation of 
Starch. Separation of Amylose and Amylopectin by Means of Some 
Nitroparaffins.”” Unpublished manuscript at the Northern Regional 
Research Laboratory. 
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